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1.  Introduction 


This  document  constitutes  the  final  report  on  a  research  study  entitled 
“Electrorheology  for  Smart  Automotive  Suspensions",  sponsored  by  the  U.S.  Army  Tank 
Automotive  Command  (TACOM)  under  contract  No.  DAAE07-92-C-R078.  The  purpose 
of  the  project  was  to  demonstrate  the  suitability  of  electrorheologicol  (ER)  technology  for 
creating  adaptive  control  of  suspension  forces  on  tracked  and  wheeled  vehicles.  Insofar 
as  ER  fluids  offer  an  inherent  electrical-to-mechanical  interface,  the  technology  is  of 
interest  wherever  the  -need  exists  for  the  active  electronic  control  of  a  mechanical 
function.  In  the  case  of  motor  vehicles,  active  suspensions  offer  a  means  of  managing  a 
broad  range  of  variables  that  otherwise  aggravate  ride  vibrations  when  only  a  fixed 
damping  device  is  installed.  Without  an  active  damping  function,  for  example,  a  fixed 
shock  absorber  must  be  designed  for  the  mean  operating  condition,  offering  too  much 
damping  when  the  vehicle  is  unladen  and,  say,  running  on  roads,  but  too  little  damping 
when  the  vehicle  is  fully  loaded  and  is  traversing  rough  terrain.  It  is  also  not  possible, 
using  a  passive  device,  to  phase  the  delivery  of  damping  forces  so  that  ride  quality  is 
optimized  even  within  a  single  cycle  of  motion.  Although  the  traditional  tools  of 
electrohydraulics  have  been  examined  for  creating  active  dampers  for  motor  vehicles,  the 
control  devices  have  typically  been  rather  complex  and  the  broad  set  of  system 
requirements  have  generally  militated  against  their  frequent  application. 

Since  ER  fluids  offer  an  inherent  medium  for  mechanical  control,  there  has  been  an 
interest  in  exploring  the  suitability  of  this  technology  for  the  active  control  of  automotive 
damping.  The  basic  concept  that  was  explored  in  this  project  involved  using  the  fluid  in  a 
device  that  directly  substituted  for  the  conventional  shock  absorber  of  the  Army’s  high- 
mobility,  multipurpose,  wheeled  vehicle  (HMMWV).  The  monitoring  of  vehicle  motions 
using  accelerometers  and  stroke  transducers  provided  the  basis  for  deriving  an  electrical 
signal,  which  was  applied  as  a  command  signal  to  the  ER-controlled  device.  The 
instantaneous  thickening  of  the  fluid  with  increasing  signal  level  produces  a  resistance  to 
motion  which  is  calculated  to  manage  the  gross  ride  vibrations  v'hich  accrue. 

The  work  reported  here  covers  a  range  of  development  tusks  including  fluid 
formulation,  development  of  an  ER  device  design,  fabrication  of  an  ER-active  damper, 
design  of  an  electrical  control  circuit,  simulation  of  ride  vibrations,  creation  of  a  quarter- 
HMMWV  test  stand,  and  laboratory  measurement  of  the  ER-damped  performance 
characteristics  over  a  range  of  physically  simulated  operating  conditions. 

This  introduction  is  followed  by  an  executive  summary  of  the  project.  The  third 
chapter  of  the  report  provides  an  overview  of  vehicle  suspension  control,  introducing  the 
notions  of  passive,  active,  and  semiactive  control,  as  well  as  the  use  of  ER  devices  as 
dampers.  Chapter  4  discusses  the  steps  that  have  been  taken  to  create  a  suitable  ER  fluid 
for  use  in  an  active  damping  device  that  is  scaled  for  the  HMMWV  application.  Chapter 
5  addresses  the  specific  control  strategy  and  the  ER  damper  design  developed  for  the 
HMMWV  in  this  project.  Chapter  C  is  a  review  of  the  physical  testing  of  the  ER  damper. 
This  includes  a  description  of  a  quarter-HMMWV  dynamometer  created  for  the  purpo.se, 
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as  well  as  a  presentation  of  test  data.  Finally,  a  simulation  study  of  the  ER  device’s 

performance  in  a  quarfer-HMMWV  model  is  reviewed  in  chapter  7. 
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2,  Executive  Summary 

The  project  combined  analytical,  design,  and  experimental  tasks  in  a  study  of 
electrorheological  (ER)  technology  applied  to  the  real-time  control  of  the  suspension  of 
the  Army’s  high-mobility,  multipurpose,  wheeled  vehicle  (HMMWV),  ER  technology  is 
based  upon  certain  fluids  which  offer  a  change  in  their  mechanical  properties  in  response 
to  a  strong  electrical  field.  The  changes  essentially  involve  resistance  to  flow  and  are 
entirely  reversible  and  very  broad  in  temporal  bandwidth.  An  ER  device,  then,  is  a 
mechanical  actuator,  damper,  valve,  or  other  element  which  contains  an  ER  fluid  and  is 
configured  to  apply  an  electric  field  across  it  for  the  sake  of  realizing  a  net  mechanical 
outcome.  In  this  project,  an  ER-based  damping  device  was  substituted  for  the  standard 
hydraulic  shock  absorber  from  an  HMMWV  suspension  to  effect  a  semiactive  damping 
control  function.  The  overall  effort  combined  the  following  elements:  (1)  the 
development  and  production  of  a  laboratory  quantity  of  ER  fluid,  (2)  characterization  of 
active  properties  of  the  ER  fluid,  (3)  development  of  an  ER  damper  suited  to  the 
mechanical  constraints  of  the  HMMWV  application,  (4)  simulation  of  the  damper  device 
as  the  active  element  in  a  quarter-HMMWV  model  (i.e.,  one  wheel  and  suspension 
assembly  supporting  one-quarter  of  the  vehicle’s  mass),  (S)  assembly  of  a  quarter-vehicle 
test  stand  for  laboratory  measurements,  (6)  fabrication  and  unit  testing  of  an  HR  damper 
and  (7)  evaluation  of  the  performance  of  the  ER  damper  by  means  of  laboratory  testing 
on  a  quarter-HMMWV  mockup. 

The  principle  behind  design  of  the  ER  damper  in  this  project  involves  the  static  stress 
capability  which  an  ER  fluid  develops  in  response  to  an  electrical  field.  That  is,  the 
electric  field  strength  causes  the  fluid  to  appear  as  a  limited-strength  solid— analogous, 
say,  to  a  stick  of  butter.  After  shearing  commences,  the  fluid’s  resistance  to  flow  depends 
increasingly  on  its  baseline  viscosity  and  less  on  the  ER  effect  of  the  static  stress 
characteristic.  Thus,  ER  devices  give  the  maximum  degree  of  control  over  a  mechanical 
outcome  if  they  can  work  the  ER  fluid  in  the  vicinity  of  the  near-static  condition.  This 
principle  was  implemented  in  the  design  of  an  HR  suspension  damper  by  incorporating  a 
screw-actuated  deyice  in  which  a  stack  of  rotary  fins  shear  the  ER  fluid  at  relatively  low 
speed,  in  resisting  the  motions  of  the  HMMW  V  suspension. 

The  prototype  ER  damper  was  sized  to  fit  within  the  coil  spring  of  the  HMMWV 
suspension,  essentially  taking  the  space  otherwise  occupied  by  the  conventional  hydraulic 
shock  absorber.  The  device  was  connected  to  a  high-voltage  amplifier  that  applies  the 
electric  field  across  alternating  pairs  of  rotary  plates,  between  which  resides  the  ER  fluid. 
Measurements  taken  on  this  device  alone,  in  a  test  stand,  showed  that  it  produced  a  static 
force  of  700  pounds  when  energized  at  a  level  that  is  seen  as  leasonable  for  long-term 
operation.  Results  obtained  by  comparing  the  ride  vibrations  achieved  with  the  ER 
damper  (as  opposed  to  the  conventional  HMMWV  shock  absorber)  showed  improved 
levels  of  isolation  with  the  ER  device.  The  primary  basis  for  improvement  stems  from 
the  fact  that  active  control  of  the  ER  damper  allows  the  precise  phasing  of  damping 
forces  so  us  to  optimally  oppose  the  motions  of  the  sprung  mass,  regardless  of  the 
immediate  state  of  the  excitation  arriving  at  the  wheels.  The  fact  that  ER  fluids  respond 
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within  a  few  milliseconds  provides  the  inhe»'ent  speed  for  accomplishing  virtually 
instantaneous  adjustments,  well  into  the  frequency  range  that  is  above  that  of  wheel-hop 
motions. 

The  conclusion  of  this  exercise  is  that  ER  technology,  while  still  unrefined  in  many 
respects,  is  an  inherently  powerful  tool  offering  completely  new  mechanisms  for  active 
control  of  automotive  functions.  The  principal  limitation  in  the  technology  at  present, 
pertains  to  the  levels  of  static  shear  stress  that  can  be  delivered  under  the  fully  energized 
condition.  Lesser  issues  pertain  to  the  settling  of  solid  particulates,  the  entrainment  of  air 
or  other  vapor  bubbles,  which  permit  electrical  arcing,  and  reduction  of  the  levels  of 
electrical  voltage  required.  The  rather  high  level  of  development  of  both  ER  fluids  and 
ER  applications  technology  that  is  proceeding  in  private  industry  suggests  that  major 
advances  in  all  of  these  areas  of  deficiency  should  be  expected. 
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3.  An  Introduction  to 
Vehicle  Suspension  Control  Techniques 

The  fundamentals  of  the  conventional,  passive  vehicle  suspension  are  reasonable  well 
understood.  Nevertheless,  good  design  requires  compromise  among  passenger  ride 
comfort,  requirement  for  suspension  working  space,  variation  in  tire/ground  contact-force 
variation,  and  vehicle  attitude  control  (Sharp  and  Crolla,  1987b). 

A  conventional  passive  suspension  consists  of  a  fixed  set  of  components,  which 
implement  predetermined  spring  and  damping  functions.  These  spring  and  damping 
functions  may  be  complex,  incorporating  discontinuous  responses  in  differing  ranges  of 
speed  and  stroke,  but  the  operation  is  entirely  mechanical  and  involves  no  external  energy 
source.  Selection  of  damper  settings  is  a  compromise  between  ride  vibration  isolation  and 
control  of  wheel-hop  resonance.  Optimal  design  can  be  made  only  under  a  defined  set  of 
conditions,  although  combinations  of  road  condition,  vehicle  loading,  tire  inflation 
pressure,  etc.  may  vary  greatly.  In  other  words,  off-optimality  is  incurred  when  the 
system  state  moves  away  from  the  assumed  design  conditions,  which  is  unavoidable. 

In  order  to  improve  the  overall  performance  of  automotive  vehicles,  suspensions 
incorporating  active  or  semiactive  components  have  been  developed.  Active  implies  that 
the  system  uses  an  external  energy  source;  semtactive,  that  adaptive  logic  is  employed 
which  changes  the  system  properties  depending  on  circumstance.  The  idea  of  using  active 
devices  to  generate  suspension  forces  dates  back  at  least  to  the  1960s,  and  this  interest  has 
remained  until  the  present,  because  of  the  large  potential  market  (Kamopp,  1990).  From  a 
control  theory  point  of  view,  the  introduction  of  active  damping  forces  would  be  expected 
to  increase  system  performance,  The  idea  that  various  active  forces  really  need  to  be 
generated  by  some  sort  of  servomechanism  is  not  particularly  attractive  in  most  cases. 
Most  force  generators  are  complex,  inefficient,  and  limited  in  frequency  response.  With 
their  necessary  power  supplies,  they  tend  to  be  much  more  expensive,  bulky,  massive, 
and  delicate  than  passive  elements.  This  led  to  the  consideration  of  semiactive  force 
generators  .such  as  the  semiactive  damper  (Alanoly  and  Sankar,  1987;  Kamopp,  1990; 
Karnopp  et  al.,  1974).  The  idea  is  to  modulate  the  dissipation  in  a  basically  passive 
damper  as  a  function  of  sensed  variables  such  as  mass  velocities. 

Passive  Suspensions 

A  conventional  passive  suspension  usually  consists  of  a  spring  and  a  shock  absorber 
(a  damper).  The  spring  can  only  store  and  return  energy  for  some  portion  of  a  suspension 
cycle,  and  the  shock  absorber  can  only  dissipate  energy.  The  damping  properties  of  the 
passive  suspension  are  constant.  There  is  no  means  of  supplying  external  energy  or  signal 
to  the  system  for  controlling  of  the  suspension  properties. 

In  order  to  cover  a  broad  set  of  conditions,  a  passive  system  can  be  designed  to  have 
parameters  varying  with  stroke  or  velocity  of  the  suspension  system.  For  example,  the 
orifice  area  varies  with  the  stroke  through  a  needle  or  plunger  with  a  variable  area.  Check 
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valves  can  be  used  to  obtain  velocity-dependent  variation,  while  pressure-relief  valves 

can  be  used  to  obtain  load-  or  velocity-dependent  properties.  However,  most  military  and 
commercial  vehicles  have  a  weight-change  factor  of  two-  or  three-to-one.  The  range  of 
the  operating  conditions  often  exceeds  the  coverage  obtained  by  varying  parameters 
passively. 

Active  Suspensions 

An  active  suspension  is  a  system  in  which  an  actuator  either  totally  replaces  the 
conventional  spring  and  damper  elements  or  acts  in  parallel  with  a  spring.  Active 
suspensions  require  external  power  to  drive  the  actuators,  which  are  normally  hydraulic 
cylinders,  to  generate  the  desired  forces  in  the  suspensions.  Such  systems  can  be  further 
classified  into  two  categories;  slow-active  and  fast-active  suspensions.  The  bandwidth  of 
a  slow-active  system  provides  for  actuation  at  frequencies  higher  than  that  of  sprung- 
mass  resonance  but  lower  than  that  of  the  fundamental  wheel-hup  vibration,  while  the 
bandwidth  of  a  fast-active  system  extends  beyond  the  domain  of  wheel-hop  frequencies. 
The  bandwidth  of  active  systems  is  mainly  liinited  by  the  type  of  actuators  used  although 
the  control  strategy  also  plays  an  important  role.  The  slow-active  systems  would  typically 
have  a  bandwidth  around  3  hz  and  would  employ  pneumatic  actuators  (Cho  and  Hedrick, 
1985),  oleo  pneumatic  actuators  (Crolla  et  ai,  1987),  d.c.  motor-driven  lead  screws 
(Sharp  and  Crolla,  1987b),  and  hydraulic  actuators  controlled  by  relatively  inexpensive 
proportional  valves  (Inagaki  et  al,  1992).  Successful  fast-active  systems  have  generally 
required  stiff  hydraulic  actuators  controlled  by  expensive,  precision  servo  valves. 

Among  the  slow-active  types  of  actuators,  pneumatic  devices  have  served  to  trim  out 
static  loads,  acting  in  parallel  with  mechanical  springs.  Such  actuators  perform  simply  as 
soft  springs  when  exercised  at  load  frequencies  that  are  beyond  their  response  bandwidth, 
such  that  the  suspension  spring  rate  stays  at  the  base  rate.  Stiff  hydraulic  actuators  as 
well  as  motor-driven  lead  screws  have  also  been  employed  in  relatively  slow-active 
systems,  being  mounted  in  series  with  mechanical  springs  (Sharp  and  Crolla,  1987b) 
when  load  trimming  is,  again,  the  objective.  Among  these  variations,  pneumatic 
actuators  are  treated  as  variable-force  actuators,  while  hydraulic  and  mechanical  lead 
screw  devices  are  seen  as  displacement  producers  changing  the  preload  on  a  conventional 
spring.  If  the  spring  is  highly  nonlinear,  as  with  most  multileaf  springs,  a  higher  preload 
also  avails  a  higher  local  value  of  spring-rate. 

The  actuators  required  for  fast-active  suspensions  are  often  not  suitable  for  use  in 
mass  produced  vehicles  such  as  cars  and  trucks  (Kamopp  and  Heess,  1991).  Since  the 
energy-input  requirements  are  high,  especially  when  addressing  high  frequency 
suspension  activity,  the  fast-active  alternative  is  adverse  to  fuel  economy  goals.  The 
Toyota  Motor  Corporation,  for  example,  recently  abandoned  their  orototype  work  on  fast- 
active  suspensions  as  impracticable  for  the  foreseeable  future  and  moved,  instead,  to 
develop  slow-active  systems  (Inagaki  et  a/.,  1992). 

In  a  slow-active  suspension,  the  active  component  generally  performs  he  leveling 
function  only,  while  its  conventional  spring-damper  subsystem  is  specially  tailored  for 
high  frequency  performance.  Therefore,  the  functions  of  attitude  control  and  isolation  are 
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separated.  Slow-active  systems  have  been  shown  to  be  competitive  with  fast-active 
systems  in  certain  constrained  applications  (Cho  and  Hedrick,  1985;  Pollard  and  Simons, 
1984;  Sharp  and  Crolla,  1987a).  However,  a  slow-active  system  may  exhibit  deficiencies 
when  the  range  of  condition  variables  is  large  and,  yet,  performance  expectations  are  high 
(Sharp  and  Crolla,  1987a).  In  particular,  slow-active  systems  fail  to  handle  wide-ranging 
conditions  when  no  mechanism  of  variable  damping  control  is  provided. 

Semiactive  suspensions 

The  most  common  semiactive  susp  nsions  are  those  whose  damping  characteristic  is 
adaptively  controlled,  although  some  have  also  entailed  direct  control  of  the  spring 
characteristic  (Mizuguchi  et  al,  1984).  In  the  spring  case,  however,  varying  the  stiffness 
without,  at  the  same  time,  effectively  varying  the  free  length  of  the  spring  produces  the 
undesirable  results  of  having  unpredictable  distribution  of  load  on  the  four  wheels  of  a 
cornering  vehicle.  Solving  this  problem  may  be  more  difficult  than  solving  those  of  a 
fully  active  system  (Sharp  and  Crolla,  1987b).  Recognizing  this  issue,  the  following 
discussion  is  limited  to  semiactive  suspensions  having  controllable  damping.  Such 
systems  can  be  fbrther  classified  into  two  sub  categories:  force-controlled  and  resistance- 
controlled  (Kamopp,  1990). 

In  the  force-controlled  semiactive  system,  the  damping  force  is  made  to  be  dependent 
only  upon  the  sign  of  the  relative  velocity;  that  is,  the  damping  force  is  rendered  more  or 
less  independent  of  the  magnitude  of  the  relative  velocity.  The  damping  force  is  raised  to 
a  step  function  in  the  tension  portion  of  the  cycle  and  set  to  zero  on  the  compression 
stroke.  This  kind  of  control  can  be  accomplished  by  (a)  using  hydraulic  shock  absorbers 
with  electromagnetically  loaded  pressure  control  valves  (Hamilton,  1985;  Kamopp  et  a/., 
1974;  Krasnicki,  1980a;  Krasnic^,  1980b),  or  (b)  incorporating  a  force  sensor  and  a  force 
control  feedback  loop  (Kamopp,  1990),  or  (c)  using  an  electromagnetic  damper 
(Kamopp,  1987). 

In  the  resistance-controlled  semiactive  system,  the  control  target  is  the  ratio  of  the 
damping  force  and  relative  velocity.  In  a  hydraulic  type  of  system,  the  damping  force  can 
be  controlled  by  changing  either  the  fluid  viscosity  or  the  orifice  area.  Control  utilizing 
fluid  viscosity  variation  will  be  effective  only  with  long-orifice  designs.  Since  the 
viscosity  of  conventional  hydraulic  fluids  can  be  controlled  only  by  thermal  regulation, 
no  feasible  control  techniques  using  conventional  fluids  exist.  Two  other  possible 
.approaches  toward  viscosity  modulation  involve  ferromagnetic  and  electrorheological 
(ER)  fluids.  The  use  of  ferromagnetic  fluids  is  generally  eliminated  from  automotive 
consideration  by  their  demand  for  high  levels  of  electric  power,  although  recent  research 
has  indicated  a  promising  new  approach  (Pinkos  et  al,  1993).  ER  fluids  offer  a  solution 
which  avoids  both  thermal  and  power-based  adjustments  in  viscous  response  and  which 
may  be  practicably  scaled  for  automotive  application. 

As  for  changing  orifice  area  to  effect  semiactive  control,  the  function  can  be 
implemented  either  continuously  or  discretely.  Thw  continuous  approach  is  to  regulate  the 
cross-sectional  area  or  opening  of  a  short  orifice  or  valve  port.  The  most  developed 
technologies  in  this  regard  involve  electrohydrauUc  servo  valves.  In  such  devices, 
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actuation  energy  is  needed  to  drive  a  valve  spool,  typically  through  either  electromagnetic 
or  piezoelectric  media.  The  discrete  approach  is  to  control  the  availability  of  flow 
passages  (each  of  which  has  a  fixed  resistance  property)  to  obtain  various  discrete 
increments  in  a  control  law  governing  total  resistance.  This  approach  can  be  realized 
using  solenoid-controlled  on-off  valves,  although  there  tend  to  be  transient  disturbances 
during  the  on-off  switching  (Kamopp  et  aL,  1974;  Margolis  and  Goshtasbpour,  1984). 

With  its  much  simpler  design,  better  efficiency  and  lower  cost,  a  semiactive 
suspension  may  compete  well  against  active  suspensions  in  practical  terms  while 
delivering  performance  that  is  much  superior  to  that  of  passive  suspensions  (Hrovat  et  ai, 
1980;  Kamopp,  1990;  Katsuda  etn/.,  1992;  Kimbrough,  1986;  Margolis,  1982;  Margolis, 
1983;  Redfield,  1991).  As  far  as  controlling  body  motion  due  to  roadway  unevenness  is 
concerned,  semiactive  systems  can  be  virtually  as  effective  as  fblly  active  systems  using 
state  variable  feedback  (Kamopp,  1990). 

Control  of  Damping  by  Means  of  ER  Fluid  Devices 

Although  a  variety  of  rather  complex  changes  in  the  mechanical  and  electrical 
properties  of  electrorheological  fluids  have  been  observed  in  response  to  an  imposed 
electrical  field,  the  primary  characteristic  upon  which  most  ER  devices  seek  to  operate  is 
the  static  shear  resistance  which  arises  exponentially  with  increasing  field  strength.  This 
phenomenon  accounts  for  the  common  observation  that  ER  fluids  can  exhibit  the 
properties  of  a  solid,  when  “energized,"  and  a  liquid  when  de-energized.  The  electrical 
Held  condition  is  imposed  upon  a  sample  of  fluid  by  containing  the  sample  between  two 
opposing  electrodes  which  are  energized  to  some  voltage.  In  order  to  realize  a 
mechanical  stress  response  from  the  fluid,  a  mechanical  shear  condition  is  established — 
typically  by  inducing  flow  of  the  fluid  between  the  electrodes  or  by  translating  one 
electrode  relative  to  the  other. 

Figure  1  shows  the  basic  characteristic  of  the  so-called  Bingham  plastic  model  of 
electrorheological  response.  As  field  strength  increases,  an  essentially  static  stress 
response  to  shear  is  observed.  Exceeding  the  static  strength  level  causes  the  fluid  to  shear 
more  or  less  according  to  the  behavior  of  Newtonian  liquids,  assuming  that  a  Newtonian 
fluid  medium  has  been  employed  as  the  liquid  phase  of  the  material. 

Since  the  imposed  field  strength  has  virtually  no  influence  on  the  shear  stress  levels 
developed  at  increasing  levels  of  shear  rate  (i.e.,  the  slopes  in  the  figure  are  virtually  the 
same  across  the  range  of  E  values,  indicating  that  viscosity,  per  se,  is  essentially 
unaffected  by  the  electric  field),  it  is  desirable  practice  to  design  devices  to  operate  in  the 
regime  of  tow  shear  rates.  By  this  approach,  variation  in  the  electrical  signal  level  more 
directly  controls  the  mechanical  (stress-based)  response,  thus  increasing  the  control 
efficiency  of  the  system  package.  The  shear  rate  regime  can  be  modified,  of  course,  by 
changing  the  nominal  gains  or  gear  ratios  which  directly  determine  either  the  flow  of  the 
fluid  through  an  electrode  gap  or  the  translational  movement  of  the  electrodes  causing 
uniform  shear  across  the  gap. 

The  particular  variety  of  semiactive  suspension  that  has  been  pursued  in  this  study 
involves  an  ER  controllable  damper  having  parallel  rotating  disks.  It  is  noted  that  a 
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Figure  li  Electrorheologic  response  per  the  classical  Bingham  model 


simple  prototype  of  ER  damper  with  concentric-cylinder  electrodes  has  also  been  tested 
by  Coulter  (Coulter,  1993)  at  Lehigh  University  within  the  same  time  frame  as  the  current 
Michigan  research. 
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4.  An  Electrorheological  Test  Fluid 

Within  this  project,  alternative  formulations  of  electrorheological  fluid  were 
examined  in  an  attempt  to  improve  upon  fluid  shear  strength  capability  and  to  reduce 
settling  of  the  particulate  component  in  ER  suspensions.  Four  alternative  fluid  categories 
were  examined,  namely,  1)  homogenous  solutions,  2)  colloidal  dispersions  or  gels,  3) 
amorphous  substituted  alumino-silicates,  and  4)  a  zeolite-based  suspension  that  was 
ultimately  used  as  the  fluid  specimen  for  testing  the  automotive  damping  device  in  this 
study.  The  consideration  of  each  will  be  discussed  below. 

Homogenous  Solutions 

ER  active  solutions  present  the  ideal  conceptual  approach  for  solving  the  problem  of 
settling  that  is  otherwise  seen  in  the  multiphase  (suspension)  materials  most  commonly 
used  as  ER  fluids.  That  is,  a  solution  is,  by  definition,  a  single  phase  system  that  is 
inherently  stable  in  terms  of  the  physical  distribution  of  its  constituents.  Recent 
discoveries  have  shown  that  solutions  of  poly  gamma  benzyl-L-glutamate  (PBLO)  and 
poly  n-hexyl  isocyanates  (PHIC)  exhibit  an  ER-active  behavior.  The  PBLOs,  although 
ER  active  in  many  solvents,  are  known  to  be  soluble  only  in  polar  solvents,  resulting  in 
relatively  high  levels  of  conductivity  and  thus  requiring  rather  high  levels  of  electrical 
power  for  sustaining  the  electric  field  condition. 

For  this  and  other  reasons,  PHIC  systems  were  examined  in  this  study,  seeking  25  to 
30  percent  concentrations  of  this  solute  in  nonpolar  solvents  such  as  xylene  and  toluene. 
The  selected  solvents  are  also  attractive  because  of  their  substantially  lower  cost  than  the 
PBLGs.  Unfortunately,  since  the  PHIC  materials  are  not  commercially  available,  it  was 
necessary  to  develop  procedures  fur  preparing  them  in  the  laboratory.  Initial  plans  were 
for  preparing  samples  oi  PHIC  from  the  monomer  form,  hexyl  isocyanate,  but  it  was 
discovered  that  only  small  quantities  were  available,  from  a  single  U.S.  vendor,  and  that 
the  highly  costly  material  could  be  obtained  only  when  ordered  with  a  nine-month  lead 
time.  Accordingly,  steps  were  taken  to  synthesize  the  monomer  in  our  laboratories  from 
more  readily  available  base  materials  (heptanoyl  chloride  and  sodium  azide).  Published 
procedures  for  this  synthesis  were  substantially  modified  to  yield  the  product  quantities 
and  the  quality  of  molecular  weight  control  that  the  ER  application  required. 
Approximately  1000  grams  of  PHIC  was  prepared,  yielding  4  liters  of  an  ER  solution. 
The  solutions  thus  prepared  were  found  to  be  substantially  ER-active  but  exhibited 
properties  which  were  poorer  than  those  of  more  common  fluids  prepared  by  other 
methods.  Since  it  was  not  possible,  within  this  project,  to  further  improve  the  materials  by 
synthetically  modifying  the  components,  the  search  for  a  suitable  ER  solution  was 
discontinued. 
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Colloidal  Dispersions/Gels 

By  a  second  approach,  steps  were  taken  to  prepare  ER  fluids  which,  under  quiescent 
or  static  conditions,  would  constitute  a  gel  which  inherently  resists  settling  of  a 
particulate  phase  of  the  materials.  The  desired  gel  formulation,  of  course,  would  exhibit 
an  ER  response  to  an  electric  field  but  would  otherwise  have  a  very  low  viscosity  when 
deformed  over  a  range  of  shear  rate  values.  Various  colloidal,  as  opposed  to  soluble, 
gelling  agents  were  used  primarily  because  the  colloidal  systems  produce  gels  that  are 
highly  shear-thinning.  By  this  approach,  the  gel  effect  would  hold  the  fluid  in  a  semi¬ 
solid  state,  at  rest,  presendng  an  essentially  zero  static  shear  strength  in  the  absence  of  an 
electric  Held.  Upon  even  very  slight  initial  loading,  the  gel  would  yield  and  thereafter 
exhibit  a  very  low  viscosity  (i.e.,  thinning)  response  to  increasing  shear  rate. 

Although  many  colloidal  agents  were  studied,  the  most  promising  was  a 
montmoiillonite  clay  modified  so  that  it  could  be  readily  dispersed  into  an  organic  liquid. 
The  ER-active  behavior  was  achieved  by  using  standard  zeolite  particles  dispersed  into 
paraffin  oil  with  1  to  2  percent  organo-clay  added.  The  fluids  had  the  appearance  of  jelly 
but  were  seen  to  be  highly  thinning,  even  at  shear  rates  down  to  0.001  per  second.  This 
approach  is  not  known  to  have  been  reported  previously.  Although  it  is  believed  to  offer  a 
promising  approach  to  achieving  stable  ER-active  suspensions,  it  was  found  that  the 
gelled  fluid  could  not  be  easily  loaded  into  the  ER  damping  device  considered  here.  That 
is,  since  the  gel-based  fluid  retains  a  static  shape,  it  was  not  possible  to  fill  a  complex 
device  with  the  fluid  without  entraining  a  great  number  of  small  air  pockets.  Entrained 
air  is  problematic  in  ER  devices  because  the  high  field  strengths  required  for  the  ER 
effect  produce  ionic  breakdown  and  arcing  across  the  bubbles.  The  arcing,  in  turn,  tends 
to  saturate  the  supply  amplifier  and  thereby  disturb  the  continuous  nature  of  the  ER 
energization  voltage.  Also,  arcing  can  lead  to  breakdown  of  many  insulating  materials 
and  to  the  introduction  of  metallic  “spatter"  which,  itself,  can  form  bridges  for  further 
short  circuit  paths  through  the  fluid.  Additionally,  optimization  of  properties  of  these 
fluids  and  techniques  for  dispersing  the  olays  into  the  liquids  are  not  straightforward  and 
will  require  substantial  research  for  their  reHnement. 

Amorphous  Substituted  Alumlno-Sillcates 

The  term,  zeolite,  applies  to  a  crystalline  fonn  of  alumino-silicates  (Al/Si)  that  were 
observed  by  University  of  Michigan  researchers  to  be  an  effective  ER-active  particulate. 
Amorphous  forms  of  Al/Si  as  the  dispersed  phase  in  ER  fluids  have  also  been  shown  to 
produce  ER  active  suspensions.  The  amorphous  materials  are  intrinsically  more 
attractive  because  the  chemistry  can  be  much  more  radically  varied  without  concern  for 
constraints  on  the  crystalline  structure  and  because  the  materials  are  much  less 
hydrophilic  than  the  crystalline  (water  entrainment  on  ER  particulates  is  undesirable 
because  of  the  higher  conductivity  level  that  accrues).  Seeking  the  combined  qualities  of 
ready  synthesis  and  dry  particles,  efforts  were  taken  to  develop  amorphous  ER  materials 
during  this  study.  While  this  task  was  carried  out  only  to  a  preliminary  level  of 
achievement,  the  effort  was  successful  in  varying  the  Al/Si  ratio  over  large  ranges  and  in 
combining  the  silicon  atoms  with  a  wide  variety  of  elements  other  than  Al.  These  include 
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boron,  germanium,  cesium,  and  silver.  Although  much  work  remains  to  be  done,  small 
fluid  samples  were  produced  showing  ER  shear  strengths  three  to  four  times  greater  than 
those  which  have  been  achieved  customarily  using  the  crystalline  materials  (zeolites). 
Research  continues  in  this  area,  given  indications  that  further  increases  in  static  stress 
capability  are  possible. 

The  ER  Fluid  Used  in  the  Damper  for  the  HMMWV 

Toward  the  conclusion  of  this  study,  a  relatively  conventional  form  of  (nominally- 
dry)  ER  fluid  was  prepared  using  commercial  zeolite  particulate  dispersed  into  a  medium 
of  transformer  oil.  The  fluid  was  quite  low  in  conductivity  and  yielded  a  shear  strength 
of  approximately  O.IS  psi  in  response  to  an  electric  field  of  3000  volts  per  millimeter. 
The  very  crucial  shear  strength  question  is  only  modestly  addressed  in  this  fluid  sample. 
Indeed  it  is  generally  accepted  that  fluid  strengths  in  excess  of  1  psi  will  be  needed  before 
ER  fluids  are  broadly  attractive  for  automotive  applications,  such  as  in  active  damping 
devices.  Nevertheless,  as  the  following  sections  of  this  document  will  show,  a  workable 
device  was  created  and  demonstrated  in  this  project. 

The  particulate  phase  of  the  test  fluid  did  tend  to  settle  when  left  on  the  shelf, 
requiring  substantial  agitation  to  redisperse  the  thick  layer  that  would  form  within  a  few 
days  or  so.  When  loaded  into  the  test  device,  the  system  was  kept  in  a  mild  state  of 
agitation  so  as  to  simply  avoid  the  settling  problem  for  the  duration  of  testing.  The 
ancillary  research  thrusts  outlined  in  the  previous  subsections  indicate  that  while  various 
approaches  exist  for  resolving  specific  shortcomings  in  ER  fluids,  a  substantial  degree  of 
further  refinement  is  still  needed. 
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5.  Design  of  a  Semiactive  Suspension  System  and  ER 
Damper  for  the  HMMWV 

In  this  chapter,  substantive  issues  of  the  ER  application  to  suspension  damping  on  the 
HMMWV  will  be  addressed.  These  issues  include  the  senUactive  control  strategy,  the 
basis  for  an  ER  damper  concept  given  certain  generic  features  of  ER  fluid  behavior,  the 
configuring  of  a  prototype  damping  device  which  is  scaled  to  the  physical  constraints  of 
the  HMMV/V  suspension,  and  simulation  of  the  dynamic  response  of  this  device  under 
conditions  that  approximate  HMMWV  operations  on  selected  roadways  from  Army 
proving  grounds.  In  the  first  subsection,  a  presentation  of  the  semiactive  control  strategy 
known  as  the  skyhook  damper  is  given.  (Later,  this  skyhook  algorithm  will  be 
incorporated  to  operate  the  ER  damper.)  Then,  a  parametric  study  of  the  basic  design 
options  establishes  the  principle  that  ER  damper  performance  is  optimized  in  devices 
which  shear  the  fluid  at  relatively  low  rates.  This  Ending  derives  from  the  observation 
that  the  electrical  field  condition  controls  a  static  shear  stress  characteristic  of  the  fluid. 
Given  this  principle,  the  discussion  proceeds  to  lay  out  a  design  in  which  the  fluid  in  a 
specialized  ER  damper  is  sheared  at  low  rate,  throughout  the  operating  range  of  the 
HMMWV,  The  device  in  question  comprises  a  screw  mechanism  which  is  advanced 
through  a  carrier  nut  which,  in  turn,  serves  to  rotate  a  number  of  disks  between  which  ER 
fluid  is  sheared.  The  screw/nut  mechanism  provides  a  m^jor  reduction  in  the  velocities 
over  which  fluid  shear  is  imposed. 

The  Semlactive  Control  Strategy 

In  this  section  the  design  of  the  semiactive  control  strategy  used  to  implement  an  ER 
damper  in  this  project,  will  be  discussed.  As  stated  earlier,  the  distinction  between  active 
and  semiactive  control  lies  in  the  field  of  force  generation.  Adjustable  damper  systems 
(such  as  the  ER-damper)  can  only  dissipate  energy  whereas  active  systems  can  also 
generate  energy.  However,  this  restriction  does  not  necessarily  imply  that  the 
performance  of  semiactive  systems  is  inferior  to  active  suspension  systems. 

Semiactive  suspension  control  is  based  on  similar  control  algorithms  as  applied  to 
active  suspension  systems.  However,  the  most  important  difference  between  active  and 
semiactive  systems  is  that  a  semiactiveiy  damped  suspension  is  only  capable  of 
dissipating  energy.  Since  the  adjustable  damper  cannot  supply  power  to  the  system,  the 
best  it  can  do  is  to  generate  no  damping  force  at  all  when  the  active  solution  requires  a 
force  that  cannot  be  generated  by  a  damper.  In  practice  this  means  that  the  ER-damper  is 
switched  to  the  lowest  setting  possible. 

In  a  broad,  simplified  sense,  improvement  in  ride  comfort  can  be  obtained  simply  by 
softening  the  suspension.  On  the  one  side,  the  suspension  spring  can  be  softened;  on  the 
other  hand  the  suspension  damper  can  be  altered  in  such  a  way  that  the  sprung  mass 
accelerations  are  reduced.  But  these  steps  must  be  taking  within  the  constraints  of 
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package  space  and  adequate  motion  control.  Semiactive  control  provides  a  good  means  to 
do  this. 

The  starting-point  of  semiactive  suspension  control  is  active  damping  based  on  an 
ideal  actuator,  situated  between  sprung  and  unsprung  mass  with  actuator  force  to  be 
specified.  The  theor>'  of  linear  quadratic  Gaussian  (LQO)  control  has  been  applied  by 
many  authors  to  determine  a  state  feedback  controller  for  active  suspension  systems.  The 
LQO  theory  is  based  on  the  n'dnimization  of  a  quadratic  performance  index.  The 
performance  of  a  quarter-car  suspension  system  can  roughly  be  assessed  quantitatively  in 
terms  of  ride  comfort  (sprang  mass  accelerations),  road  holding  (variations  of  the  vertical 
tire  load)  and  suspension  working  space.  These  three  criteria  can  be  used  to  setup  a 
quadratic  performance' criterion.  The  relative  importance  of  the  different  performance 
parameters  can  be  expressed  by  weighting  factors.  After  having  chosen  the  weighting 
factors,  the  optimal  control  algorithm  will  select  a  state  feedback  gain  matrix  such  that 
the  quadratic  cost  function  is  minimized.  This  feedback  can  either  be  a  full  state  feedback 
or  a  partial  state  feedback.  The  full  state  feedback  considers  all  the  states  of  the  system 
(such  as  displacements  and  velocities  of  sprang  and  unsprung  masses)  and  uses  this 
information  to  compose  an  active  suspension  force  applied  between  axle  and  vehicle 
body.  On  the  other  hand  the  limited  state  feedback  looks  only  at  particular  states  (such  as 
the  sprung  mass  velocity).  The  limited  state  feedback  is  favorable  because  fewer  states 
have  to  be  measured  when  the  controller  is  implemented  in  a  real  vehicle.  Venhovens 
(1993)  has  shown  that  the  performance  of  limited  state  feedback  systems  in  terms  of  gain 
in  ride  comfort  is  almost  equal  for  both  types  of  feedback. 

The  ultimate  way  to  reduce  sprung  mass  accelerations  is  to  apply  an  actuator  force 

which  is  merely  proportional  to  the  sprang  mass  absolute  velocity  according  to 

(1) 

where  represents  the  so-called  skyhook  damping  constant.  Since  only  the  sprung  mass 
absolute  vertical  velocity  is  used  in  equation  (1),  thi3  type  of  feedback  is  an  example  of 
limited  state  feedback.  The  active  skyhook  configuration,  together  with  a  real  skyhook 
system  and  a  semiactive  damper  system  with  skyhook  control  are  shown  in  figure  2. 


m, 


/////// 

F iyn 

Figure  2.  Skyhook  control,  from  active  to  semiactive 
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Figure  3.  Frequency  response  ftincttons  of  three  suspension  systems 
(a)  tire  deflection  (b)  sprung  mass  acceleration 

(c)  suspension  working  space  (d)  sprung  mass  displacement 
(e)  actuator/damper  force 

Although  there  is  a  substantial  difference  between  a  real  skyhook  system  and  the 
active  counterpart  (the  actively  generated  force  acts  on  both  sprung  and  unsprung  mass; 
the  real  skyhook  damper  force  acts  only  on  the  sprung  mass),  the  frequency  response 
functions  of  both  systems  show  great  correspondences.  The  five  response  functions  of 
interest  are  plotted  in  figure  3. 

The  modulus  of  the  tire-deflection  transfer  function  (figure  3)  of  the  active  system 
with  skyhook  control  differs  from  the  curve  valid  for  a  real  skyhook  system.  The  reason 
is  the  point  of  action  of  the  damper/actuator  force.  Not  totally  unexpected  is  the  presence 
of  a  huge  resonance  peak  nearby  12.6  hz  (unsprung  mass  natural  frequency).  Since  the 
motion  of  the  unsprung  mass  is  not  damped  properly  in  case  of  the  skyhook  damper  or 
active  system  with  skyhook  control,  a  resonance  peak  is  the  result, 

A  small  but  not  insignificant  difference  is  the  height  of  the  resonance  peak.  In  case  of 
the  tmc  skyhook  system  the  wheel-hop  mode  is  nut  damped  at  all,  resulting  in  an  infinite 
modulus  at  that  specific  frequency.  The  active  system  with  skyhook  control  however  has 
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a  bounded  finite  amplitude.  Although  small,  the  sprung  mass  velocity  is  not  equal  to  zero 
near  the  wheel-hop  natural  frequency.  This,  combined  with  the  fact  that  the  actuator  force 
depends  solely  on  the  sprung  mass  velocity,  and  that  this  force  acts  on  both  sprung  and 
unsprung  masses,  results  in  a  finite  modulus. 

With  exception  of  the  area  near  the  unsprung-mass  natural  frequency,  both  skyhook 
versions  (true  one  and  active  approximation)  offer  a  large  improvement  in  vibration 
isolation  (graph  b)  over  the  entire  frequency  range  compared  with  the  passive  system 
with  a  fixed  suspension  damping.  Anticipating  the  test  result  of  the  quarter-car  model 
with  ER-damper,  an  improvement  in  ride  comfort  and  low  frequent  tire  load  variations  is 
to  be  expected  at  the  cost  of  a  degradation  of  tire  load  variations  in  the  range  of  the 
wheel-hop  natural  frequency. 

Generic  Design  Options  and  Parameters  for  the  ER  Damper 

In  approaching  the  design  of  an  ER-based,  semiactive  suspension  control,  the  ER 
component  can  be  built  as  either  a  flow-mode  (flgure  4a),  mixed-mode  (figure  4b),  or 
shear-mode  (figure  4c)  type  of  damper.  The  configurations  illustrated  in  figure  4  are 
essentially  generic  concepts  that  contain  simplified  features  to  facilitate  the  design 
analysis.  In  practical  hardware,  for  example,  the  cylinder  is  likely  to  be  single-ended, 
rather  than  double-ended,  and  the  ER  control  valve  or  area  is  likely  to  be  comprised  of 
multiple  parallel  or  concentric  plates.  Among  the  three  conceptual  modes,  the  flow-mode 
damper  is  most  similar  to  the  traditional  shock  absorber  except  that  it  replaces  the 
conventional  orifice  with  an  ER  control  gap  or  valve.  The  source  of  the  damping  force  in 
this  mode  is  exclusively  from  flow-induced  pressure  drop  across  the  piston.  The  ER  valve 
can  be  placed  outside  the  cylinder  (as  in  figure  4a)  or  within  the  piston  component.  For 
the  sake  of  compact  packaging,  the  external  approach  could  be  implemented  by  forming  a 
flow  chamber  along  the  outside  of  the  cylinder  wall  (Petek,  1992). 

In  the  mixed-mode  damper,  the  ER  control  is  realized  in  the  gap  between  the  piston 
side-wall  and  the  cylinder  liner.  The  fluid  shear  stress  at  the  shear  surface  also  contributes 
to  the  damping  force,  while  the  flow-induced  pressure  drop  is  still  present  as  a  force 
component  on  the  piston. 

If  a  large  by-pass  port  is  cut  through  the  piston  of  a  mixed-mode  damper,  the  device 
becomes  a  simple  shear-mode  element.  Since  the  fluid  can  flow  freely  through  the  by¬ 
pass,  the  pressure  drop  across  the  piston  becomes  negligible  and  only  shear  stre.sses 
developed  along  the  wall  contribute  to  the  damping  force. 

Our  analysis  (see  appendix  A)  indicates  that  the  shear-mode  damper  achieves  better 
dynamic  response  and  more  effective  control  than  the  flow-  or  mixed-mode  types. 
Because  certain  principles  in  ER  design  support  the  prospect  of  a  package  that  is 
mechanically  simple,  durable,  low  in  cost,  and  exhibits  a  broad  bandwidth,  it  appears  that 
an  ER  semiactive  suspension  may  be  superior  to  conventional  types  of  semiactive  and 
active  suspensions  for  automotive  applications.  Although  the  fluid’s  behavior  is  complex, 
it  has  also  been  sufficiently  well  defined  at  this  point  that  its  codification  on  cheap  and 
reliable  memory  devices  is  straightforward.  Improvement  in  mechanical  strength  of  the 
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(c)  Shear  -mode  (d>  ER  control  gap 


Figure  4.  The  generic  concepts  of  the  three  modes  of  electrorheological  dampers 

fluid’s  energized,  shear-stress  response  remains  the  most  important  factor  in  determining 
whether  or  not  ER  devices  will  become  practical  in  this  application. 

Selected  Design  of  the  ER  Shock  Absorber 

A  shear-mode  damper,  as  indicated  in  the  above  section,  achieves  better  dynamic 
response  and  more  effective  control  than  the  flow-  or  mixed-mode  types  and  is  thus 
selected  as  the  damper  mode  for  this  project.  The  general  concept  of  the  proposed  design 
is  shown  in  figure  S.  The  control  in  the  ER  damper  is  realized  through  a  medium  of  an 
ER  fluid.  The  shear-mode  damper  uses,  instead  of  translational  plates,  a  plurality  of 
rotational  shearing  disks  (hot  electrodes  and  cold  electrodes),  which  form  a  total  control 
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(or  shearing)  surface  area  that  is  constant  and  independent  of  the  relative  position  of  the 
actuator,  resulting  in  a  good  control  and  a  compact  size.  A  translation-to-rotation  device 
converts  the  translational  motion  of  the  vehicle  vibration  into  the  rotational  motion  of  the 
shear-mode  damper,  and  vice  versa. 

The  actual  ER  damper  developed  in  this  project  is  show  in  figure  6.  The  specific 
design  of  this  device  is  illustrated  in  figure  7.  In  this  design,  thirty-nine  ground 
electrodes,  separated  from  each  other  by  spacers,  are  mechanically  fixed  and  electrically 
connected  to  a  stator  through  stator  feather  keys.  The  stator  is  grounded  electrically  and 
fixed  rigidly  onto  the  unsprung  mass.  Forty  elevated-voltage,  or  “hot,"  electrodes,  also 
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Figure  <i.  The  ER  damper  devdopcd  fm*  the 
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Figure  7.  Assembly  of  the  ER  damper 
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separated  from  each  other  by  spacers,  are  mechanically  attached  to  a  rotator  through 
rotator  feather  keys,  and  electrically  connected  to  the  high  voltage  through  a  spring- 
loaded  carbon  brush.  The  rotator  is  mechanically  supported  by  top  and  bottom  bearings 
that  can  sustain  both  radial  and  axial  loads.  The  rotator  is  made  from  glass  epoxy  0-1 1,  a 
dielectric  material  with  high  mechanical  strength.  Interelectrode  gaps  are  filled  with  an 
ER  fluid  through  top  and  bottom  ports. 

The  translational  relative  motion  is  converted  through  a  ball  screw  mechanism,  which 
is  tightly  connected  on  top  of  the  rotator  with  a  joining  nut.  The  use  of  the  ball  screw 
renders  smooth  power  and  motion  transmission.  The  elevated  voltage  can  be  supplied  to 
the  rotator  through  a  carbon  brush  at  the  bottom  instead  of  the  fixed  connection  at  the  top 
as  shown  in  figure  5.  The  complete  set  of  drawings  of  the  ER  damper  is  contained  in 
appendix  B. 

Target  Parameters  for  ER  Shock  of  TACOM 

A  conventional,  rear  suspension  system  of  the  HMMWV  is  illustrated  in  figure  8.  A 
shock  absorber,  made  by  the  Monroe  Auto  Equipment  Company,  is  mounted  inside  a  coil 
spring.  To  make  minimum  changes  to  the  existing  suspension  system,  the  ER  damper  is 
also  situated  inside  the  mechanical  spring,  and  its  outer  diameter  is  therefore  limited  by 
the  inner  diameter  of  the  spring,  which  varies  between  1 1.5  cm  (near  the  bottom  seat  of 
the  spring)  and  12.6  cm  along  the  axis.  The  selection  of  the  outer  diameter  of  the  ER 
damper  (equal  to  the  outer  diameter  of  the  stator)  is  further  influenced  by  the  spring 
budding  and  some  misalignment  between  the  spring  and  shock  center  lines.  The  outer 
diameter  of  the  Monroe  shock  is  6.985  cm  (2.75  in).  For  a  quarter-HMMWV,  the  sprung 
and  unsprung  masses  are  803  kg  (1769  lbs)  and  98  kg  (217  lbs),  respectively.  The 
suspension  spring  stiffness  measured  at  the  wheel  spindle  is  63,528  N/m  (4,356  Ibs/ft  or 


Figure  8.  The  conventional  rear  suspension  system  of  an  HMMWV 
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363  Ibs/in).  The  tire  spring  stiffness  is  204,394  N/m  (14,015  Ibs/ft).  The  vehicle  heave 
natural  frequency  is  1.3  hz,  while  the  wheelhop  frequency  is  7.9  hz.  The  vehicle  pitch 
natural  frequency  is  1.8  hz,  which  is  ignored  here  because  only  a  quarter  car  model  is 
considered. 

The  metal-to-metal  length  of  the  Monroe  shock  is  32.72  -  32.1 1  to  42.16  -  41.55  cm 
(12.88  -  12.64  to  16.60  -16.36  in),  and  the  metal-to-metal  stroke  is  9.4488  cm  (3.72  in). 
The  overall  stroke  is  5.3  cm  (2.1  in).  The  length  of  the  body  of  the  shock  absorber  is 
about  21.3  cm  (8.38  in).  Because  of  the  geometric  layout,  the  displacement  and  force  at 
the  wheel  centerplane  (Xwheel  and  Fwheel)  w®  related  to  those  at  the  shock  absorber 
(Xghock  and  Fshock)  approximately  in  the  following  way: 

^shock  ®  0*48  Xwheel  (2) 

Fshock  ■  2.08  Fwheel 

A  conventional  shock  has  bilinear  characteristics.  It  offers  more  damping  during 
extension  than  during  compression.  Also,  the  shock  has  a  pressure  relief  feature  which 
limits  the  damping  force  to  a  maximum  value  when  the  relative  velocity  creates  a  certain 
relief  pressure.  Some  of  the  base  data  are  presented  in  table  1.  Relief  velocities  and  limit 
damping  forces  at  the  wheel  spindle  were  taken  directly  from  the  vehicle  test  data  sheet 
from  TACOM.  Using  the  above  two  equations,  the  spindle  velocities  and  forces  were 
extrapolated  to  the  shock  velocities  and  forces,  which  do  not  always  agree  with  the  shock 
data  provided  by  Monroe.  Also  damping  forces  are  far  from  linear  within  the  limit 
velocities  and,  thus,  their  corresponding  damping  coefficients  are  not  constant.  Their 
simple  arithmetic  average  values  were  calculated.  Damping  ratios  were  also  calculated, 
using  the  sprung  mass,  the  spring  constant,  and  the  average  damping  coefficients.  As 
expected,  the  damping  ratio  for  extension  ( =  0.74)  is  higher  than  the  damping  ratio  for 
compression  ( s  0.25). 

Based  on  the  above  datci,  we  chose  15,000  N  (3,374  Ibf)  as  the  target  design  damping 
force  in  the  ER  damper. 

Number  of  Electrode  Plates 

The  influence  of  the  inertia  force  is  small  in  a  shear-mode  damper,  such  that  the 
damping  force  in  the  strut  can  be  approximated  with  a  Newtonian  damping  force  (Fjn) 
and  an  ER  damping  force  (Fjer),  corresponding  to  the  effects  of  the  plastic  viscosity  (T|) 
and  the  yield  shear  stress  (xy),  respectively,  in  the  Bingham  model.  The  two  damping 
forces  are  derived  us 

(4) 

ehD^an  k 

and 

FaER^?52i!i^Ty  (5) 

l^''l  3cDptan  X, 


with 
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Table  1.  Basic  parameters  of  the  HMMWV  shock  absorbers 


Compression 


At  the  wheel  At  the  shock 


0.76  m/s 
2.48  ft/s 


Relief  velocity 


Limit  damping  2,667  N 

force  at  and  ^  Ibf 

beyond  relief 
velocity 


Damping  2,348  to  4,696 
coefficient  range  '  N  s/m 
161  to  322 
Ibf  s/ft 


Average  damping 
coefficient 


Average  damping 


Monroe  shock  limit 
damping  force  and 
velocity 


3,522  N  s/m 
242  Ibf  s/ft 


0.36  m/s 
.19  ft/s 


5,547  N 
1,248  Ibf 


1,190-714  Ibf 
at  1.11  ft/s 


Extension 


At  the  wheel  At  the  shock 


0.76  m/s 
2.48  ft/s 


10,209  to 
10,923  N  s/m 
700  to  749 
Ibf  s/ft 


10,566  Ns/m 
725  Ibf  s/ft 


0.36  m/s 
1.19  ft/s 


16,643  N 
3,744  Ibf 


2,608-1,564 
Mat  1.11  ft/s 


(6) 

dt  dt 

cosA  •  ^ 

- lanAt 

(7) 

ftanX>f  cositi 

where  n  is  the  number  of  hot  electrodes,  /»  (  *  0.001  m)  the  gap  size  between  the 
electrodes,  t  the  time,  Av  the  strut  deflection  velocity,  the  sprung  mass  displacement, 
xy  the  unsprung  muss  displacement,  ro  (»  0.0483  m)  and  ri  (=  0.0285  m)  the  outer  and 
inner  radii  of  the  effective  shearing  surface  of  the  electrodes,  respectively.  The  parameter 
Dp  (  =  0.0254  m)  is  the  pitch  diameter,  (|>  (standard  at  29  degrees)  the  thread  angle  ,  X  (  = 
17.657  degrees)  the  lead  angle,  e  (  =  90  percent)  the  ball  screw  efficiency,  and  /  (  = 
.033)  the  friction  coefficient  for  the  screw-nut  mechanism.  The  variables  vs  and  vy  are 
the  sprung  and  unsprung  mass  velocities,  respectively,  and  Av  is  the  strut  deflection 
velocity.  The  Newtonian  and  EK  damping  forces  here  are  completely  uncoupled  and 
additive.  While  this  approach  provides  only  a  first-order  estimation,  it  is  convenient  and 
insightful  for  a  preliminary  parametric  study.  However,  a  complete  solution  of  the 
problem  will  demand  a  full  solution  of  the  fluid  flow  problem  (Lou  et  ai,  1993). 

If  FriER  alone  is  to  reach  15,000  N  and  if  Xy  =  6,895  Pa  (1  psi),  n  must  be  21  according 
to  equation  (6).  If  xy  =  0.5  psi,  then,  n  =  42.  In  the  current  design,  40  hot  electrodes  are 
chosen. 

Damping  torque  on  the  rotator  Tj  is  related  to  the  damping  force  (Fj)  as  follows 
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Tjj  =  e^an  XFd 

m 

where  Fd  *  15.000  N,  Td  *  54.57  N-m. 


(8) 


Material  of  the  Rotator 


The  material  used  for  the  rotator  must  satisfy  both  mechanical  and  electrical 
constraints.  The  maximum  shear  stress  is 


tKD-'  -  d*) 


(9) 


where  D  and  d  are  the  outer  and  inner  diameters  of  the  rotator,  respectively.  In  this 
design,  D  «  39.7  mm  and  d  s  27.99  mm.  Therefore, 

Sv«5.9xlO<^N/m2  ^O) 

1'he  peak  intensity  of  the  electric  field  requires  that  the  rotator  be  able  to  withstand  an 
8  kv/mm  field  condition.  To  satisfy  the  combined  requirements,  a  glass  epoxy  0-11  was 
selected.  This  material  provides  a  shear  strength  of  1.517  x  10^  N/m^  (22,000  psi)  and  an 
electric  strength  of  35.4  kV/mm  (900  V/mil). 


Rotator  Feather  Key 

According  to  moment  balance, 

Td-ricOtdtknicn  (11) 

where  rih  is  inner  radius  of  the  hot  electrode  plate  (■  1.877  in  ■  0.0477  m);  r,  shear 
stress,  tj,  the  hot  electrode  plate  .thickness  («  0.001  m);  ,  the  depth  of  the  rotator  key  {» 

0.197  in  ■  0.005);  the  number  of  the  keys  (b  2).  The  maximum  imposed  torque 
results  in  a  shear  stress  on  the  key  of: 

a  Bs  Td/ri(t<it|cnkn  *  2.86  x  1 0®  Wrrfi  / 1  o') 


Stator  Feather  Key 

Each  stator  key  is  fixed  by  eight  screws  for  easy  manufacturing.  Each  screw  will  bear 
a  shear  stresses  according  to  the  relationship: 

TdB8r,a,n,7id?/4  (13) 

where  rg  is  the  outer  radius  of  the  ground  electrode  plate  (s  4.13  in  s  0.105  m);  n^,  the 
number  of  the  stator  keys  (s  2);  dg,  the  screw  diameter  (»  0.1 12  in  »  0.00284  m).  The 
maximum  shear  .stress  on  each  stator  key,  Oj,  is  calculated  to  be  5.39  x  lO''  N/m^,  which 
is  smaller  than  the  allowed  stress  of  35  x  10**  N/m^. 
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Filling  the  Damper  with  £R  Fluid 

In  order  to  avoid  arcing  when  high  voltages  are  applied  across  the  plates  of  the 
damper,  it  is  important  to  insure  that  little  or  no  air  is  entrained  in  the  ER  fluid. 
Therefore,  a  special  effort  was  undertaken  to  fill  the  damper  with  fluid. 

The  fluid  sample  was  first  put  into  a  closed  vessel  and  subjected  to  a  29-in  hg. 
vacuum  for  a  period  of  twenty-four  hours.  Still  under  a  vacuum,  the  fluid  was  agitated  to 
insure  an  even  distribution  of  the  particulate  material,  The  fluid  was  then  immediately 
forced  into  the  damper  under  three  to  five  psi  pressure.  In  this  process,  the  fluid  was 
withdrawn  from  the  bottom  of  the  external  vessel  and  fed  into  the  damper  through  its 
bottom  port  so  that  air  was  not  drawn  into  the  system.  Air  was  simultaneously  evacuated 
from  the  damper  through  the  top  port. 

In  all  of  the  subsequent  operations  of  the  damper,  the  internal  fluid  space  of  the 
damper  was  plumbed  to  a  small  bladder-accumulator  containing  the  same  ER  fluid.  By 
this  means,  ER  fluid  was  maintained  at  about  three  psi  gage  pressure.  Thus,  the  fluid  in 
the  damper  was  maintained  at  a  positive  pressure,  preventing  air  ingestion  through  the 
seals  of  the  damper  while  permitting  expansion  and  contraction  with  heating  and  cooling. 

Whenever  the  damper  was  not  active,  the  sealed  system  of  damper  and  external 
accumulator  was  removed  from  the  testing  apparatus  and  placed  in  a  tumbler  in  order  to 
prevent  settling  of  the  particulates  in  the  ER  fluid.  (The  tumbler  was  simply  a  tube  lined 
with  foam  padding,  turning  slowly  in  a  lathe.) 


27 


6.  Testing  the  ER  Fluid  Damper 

Three  types  of  tests  were  conducted  on  the  semiactive  ER  damper  developed  by 
UMTRI.  The  initial  tests  involved  a  very  simple  series  of  static  measurements  to 
quantify  the  force-to-voliage  behavior  of  the  new  damper.  These  tests  were  followed  by 
two  series  of  more  sophisticated  dynamics  tests.  In  one,  the  damper  was  tested  alone  in 
order  to  verify  its  basic  properties.  The  third,  more  elaborate  set  of  tests  was  conducted 
with  the  damper  installed  in  a  quarter-HMMWV  dynamic  test  facility.  These  tests 
measured  the  ride  isolation  properties  of  the  HMMWV  suspension  using  both  the 
semiactive  ER  damper  and  the  conventional  damper  normally  used  in  the  HMMWV 
suspension. 

This  chapter  is  divided  into  four  sections  which  cover  (1)  the  static  test  series,  (2)  the 
dynamic  test  facility,  (3)  dynamic  tests  of  the  damper  alone,  and  (4)  the  quarter- 
HMMWV  series  of  tests. 


Static  Tests  of  the  ER  Damper 

The  static  tests  were  done  simply  to  measure  the  relationship  between  applied  control 
voltage  and  the  static  force  developed  by  the  damper.  As  noted  previously,  the  primary 
quality  of  the  shear-mode  type  of  damper  is  that  it  produces  a  force  in  response  to  control 
voltage  (i.e.,  the  Bingham  shear  response)  but  is  rather  insensitive  to  stroke  velocity  (i.e., 
the  viscous  shear  response).  It  follows  that  a  static  characterization  of  the  control 
voltage-to-force  relationship  of  the  device  can  serve  as  a  primary  description  of  the 
behavior  of  the  device. 

The  static  tests  were  very  simple.  The  body  of  the  damper  was  fixed  on  the  test 
bench  and  a  control  voltage  was  applied.  A  torque  wrench  was  then  used  to  rotate  the 
ball-screw  nut.  The  measured  torque  was  converted  to  an  equivalent  force  for  the  ball- 
screw  assuming  one  hundred  percent  mechanical  efficiency.  The  process  was  repeated  at 
several  voltages. 

The  results  of  two  repeats  of  this  simple  test  are  shown  in  figure  9.  Two 
characteristics  of  this  plot  are  significant.  Firstly,  the  damper  exhibits  a  coulomb  friction 
response  of  about  50  pounds  force  with  no  voltage  applied.  This  force  derives 
exclusively  from  the  rotational  friction  of  the  fluid  seals  amplified  by  the  mechanical  gain 
of  the  ball-screw  hardware.  Secondly,  the  force  response  of  the  system  is  approximately 
proportional  to  the  2.5  power  of  the  applied  voltage.  This  is  similar  to  previous 
observations  of  the  behavior  of  ER  fluid  systems  (Lou,  1990). 

The  data  of  figure  9  have  been  fitted  with  two  curves.  One  curve  is  a  quite  good  fit 
which  reflects  the  coulomb  friction  quality  of  the  system  as  well  as  the  voltage  response. 
The  second  is  a  poorer  fit  in  that  it  ignores  the  friction  response  and  approximates  only 
the  voltage  response.  (That  is,  this  curve  passes  through  the  origin  of  the  plot.)  This 
latter  function  was  used  us  a  basis  for  the  semiactive  control  algorithm  in  the  quarter- 
HMMWV  test  series. 
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Figure  9.  The  static  response  of  the  ER  damper 


The  Dynamic  Test  Facility 

In  order  to  conduct  in  situ  dynamic  tests  of  the  ER  damper  us  a  semiactive  shock 
absorber,  UMTRI  constructed  a  quarter 'HMMWV  dynamic  road  simulator  facility. 
TACOM  provided  the  necessary  HMMWV  suspension  and  frame  components.  These, 
along  with  the  necessary  supports,  instrumentation,  and  electrohydraulic  servo  road 
surface  exciter,  made  up  the  quarter-HMMWV  road  simulator.  The  facility  was  used  to 
measure  the  ride  response  of  the  quarter-HMMWV  to  road  profiles  provided  by  TACOM. 
Test  were  conducted  with  the  standard  damper  used  on  the  HMMWV  as  well  as  with  the 
ER  damper. 

In  addition  to  using  this  equipment  for  fuil  scale  quarter-HMMWV  tests,  a  subset  of 
the  servo  actuator  and  instrument  systems  were  used  for  a  small  series  of  dynamic  tests 
on  the  ER  damper  alone.  The  servo  actuator  was  used  to  stroke  the  damper  in  a 
controlled  manner.  Tests  were  conducted  to  determine  the  response  of  the  damper  to 
stroke  velocity  and  to  step  changes  in  control  voltage. 

The  remainder  of  this  section  will  describe  the  dynamic  test  facility.  The  dynamic 
tests  themselves  will  be  discussed  in  sections  to  follow. 
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The  Mechanical  Facility 

The  quarter-HMMWV  road  simulator  facility  is  pictured  in  figure  10.  Figure  1 1  is  a 
sketch  of  the  facility  in  which  mtniy  of  the  items  mentioned  in  the  following  discussion 
are  labeled  for  the  reader’s  convenience. 

The  facility  is  based  on  a  small  bed  plate.  Two  pairs  of  channel-section  uprights  are 
mounted  on  the  plate.  The  shorter  pair,  immediately  fore  and  aft  of  the  HMMWV  tire, 
support  the  vertical  linear  bearings,  which  guide  the  motion  of  the  simulated  road  surface. 
The  taller  pair  of  uprights  support  the  vertical  guide  bearings  for  the  motion  of  the  vehicle 
sprung  mass. 

This  “vehicle”  is  one  rear  quarter  of  an  HMMWV,  The  tire,  wheel,  hub,  and 
suspension  members  (except  for  the  modified  shock  bracket)  are  all  standard  HMMWV 
parts.  The  suspension  is  attached  to  a  shortened  HMMWV  frame  rail  and  standard  frame 
cross  members.  (All  the  standard  HMMWV  parts  in  figure  10  are  painted  black.)  A  gray 
colored  ballast  frame  is  mounted  on  the  standard  frame  members.  This  frame  supports 
three  500-pound  iron  blocks,  bringing  the  total  mass  of  the  quarter-HMMWV  to  2,490 
pounds,  thus  simulating  a  fully  loaded  HMMWV.  Note  that  the  ballast  iron  is  mounted 
directly  above  the  tire.  This  helps  locate  the  center  of  gravity  of  the  simulated  vehicle 
closer  to  the  wheel  plane,  greatly  reducing  the  overturning  moment  applied  to  the  linear 
bearings,  which  provide  the  vertical  guidance  for  the  sprung  mass. 

The  dynamic  road  simulator  is  located  beneath  the  HMMWV  tire.  The  road  surface 
plate  is  mounted  on  linear  bearings,  which  accommodate  the  lateral  motion  of  the 
HMMWV  tire.  (Because  of  compromises  forced  by  ground  clearance  requirements,  the 
HMMWV  suspension  geometry  produces  an  unusual  amount  of  lateral  motion  of  the  tire 
when  the  suspension  deflects  vertically.)  These  bearings  sit  on  an  intermediate  plate, 
which  is  guided  by  the  vertical  linear  bearings  mentioned  previously. 

Vertical  support  for  the  road  is  provided  by  three  cylinders.  The  outer  two  are  air 
cylinders,  which  provide  support  for  the  static  load  of  the  system.  The  center  cylinder  is 
an  hydraulic  servo  actuator,  which  provides  the  dynamic  component  of  the  road  force. 
All  three  cylinders  are  mounted  atop  an  air  manifold. 

The  air  cylinders  are  plumbed  to  the  air  manifold  through  a  large  (2-in  diameter)  hole 
bored  directly  down  through  the  cylinder  tail  block  and  on  into  the  air  manifold.  In  turn, 
the  manifold  is  plumbed  through  a  short  2-in  pipe  to  a  large  air  tank.  The  tank  is  supplied 
with  pressure-regulated  air.  This  large  volume,  low-flow-resistance  system  results  in  the 
air  cylinders  providing  a  virtually  constant  force  to  the  road  even  during  large  and  rapid 
road  motions.  The  air  pressure  is  adjusted  so  that  this  static  force  is  approximately  equal 
to  the  weight  of  the  quarter-HMMWV. 

Because  of  the  air  system  support,  the  hydraulic  actuator  need  supply  only  the 
dynamic  portion  of  the  road  forces.  Relieving  the  servo  of  the  static  support  substantially 
reduces  the  power  requirements  of  the  servo  system,  allowing  higher  performance  with 
smaller,  less  expensive  components.  The  2-1/2-in  diameter  hydraulic  cylinder  was 
supplied  through  a  pair  of  parallel  plumbed  lO-gpm  MTS  2S2-series  servo  valves. 
(These  are  actually  modified  Moog  76-scries  valves.)  Hydraulic  power  was  provided  by 
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Figure  10.  Photograph  of  the  dynamic  test  facility 
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a  10-gpm,  3000  psi  supply  augmented  with  a  1-1/4-gaIlon  supply  accumulator.  These 
components  formed  a  system  whose  small  displacement  response  was  flat  to 
approximately  40  hz.  A  maximum  velocity  of  20  in/sec  was  attainable  within  the 
accumulator  reserve.  Sustainable  mean  velocity,  established  by  pump  rate ,  was  about  18 
in/sec. 

Instrumentation  and  Control 

Figure  12  presents  a  block  diagram  of  the  instrumentation  and  control  electronics  of 
the  quarter-HMMWV  road  simulator.  The  system  provides  for  the  three  basic  functions 
of  (1)  ER  damper  control,  (2)  data  acquisition,  and  (3)  road  actuator  control.  The  heart  of 
the  system  was  a  Macintosh  Quadra  800  computer  coupled  to  UMTRI's  Engineering 
Research  Division  (ERD)  data  acquisition  system.  This  basic  system  was  augmented 
with  a  high>voltage  amplifier  for  generating  the  kilovolt'level  signal  required  for  damper 
control,  high'Voltage  signal  isolators  for  protecting  the  computer  circuits  in  case  of 
arcing,  a  servo  amplifier,  thermocouple  ampllHer,  and  an  array  of  appropriate  sensors. 
Details  of  these  elements,  in  the  context  of  the  three  basic  functions,  are  presented  in  the 
following  three  subsections. 

ER  Fluid  Damper  Control 

The  computer  controller  of  the  quarter-HMMWV  simulator  determined  the  electrical 
voltage  sent  to  the  ER  damper.  This  voltage  command  (Vq)  was  recalculated  and  updated 
100  times  per  second.  The  calculated  command  was  output  by  the  computer  through  a 
high-voltage  isolator  into  a  high-voltage  amplifier.  (See  the  next  section  for  a  discussion 
of  these  two  elements.)  The  actual  voltage  signal  (v)  —  up  to  3000  volts  —  was  supplied 
to  the  ER  damper  by  this  amplifier. 

In  concept,  any  of  the  forces  or  motions  transduced  by  the  many  instruments  of  the 
simulator  could  be  used  by  the  computer  to  calculate  the  desired  ER  damper  voltage.  In 
these  tests,  only  two  velocity  signals  were  actually  used.  These  were  the  absolute  vertical 
velocity  of  the  sprung  mass  (Vs)  and  the  relative  velocity  of  the  sprung  and  unsprung 
masses  (V(s.u)),  that  is,  the  damper  stroke  velocity.  These  velocities  were  used  in  two 
control  algorithms,  one  which  caused  the  ER  damper  to  behave  as  a  skyhook  damper  and 
one  which  caused  it  to  behave  as  a  conventional  shock  absorber. 

The  skyhook  damper  was  discussed  in  chapter  5  of  this  document.  To  briefly  review: 
The  concept  is  one  of  a  linear  damper  connected  between  the  sprung  mass  of  the  vehicle 
and  an  imaginary  hook  in  the  sky  whose  vertical  velocity  is  defined  to  be  zero.  Thus,  the 
damping  forces  are  proportional  only  to  the  absolute  vertical  velocity  of  the  sprung  mass, 
not  to  the  relative  velocity  of  the  sprung  and  unsprung  masses  as  is  the  case  with  a 
conventional  damper. 

In  practice,  the  real  damper,  which  produces  these  idealized  forces,  is  actually 
connected  between  the  sprung  and  unsprung  masses.  Since  the  damper  is  passive  (that  is, 
it  can  produce  forces  only  in  opposition  to  stroke  velocity),  there  are  time  periods  when  it 
is  not  possible  for  the  damper  to  generate  the  desired  damping  forces.  For  example,  when 
1  he  sprung  mass  is  moving  upward,  the  desired  force  is  downward,  but  the  damper  can 
i  pply  a  force  downward  on  the  sprung  mass  only  when  the  relative  velocity  of  the  sprung 
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Block  diagram  of  the  dectrMiic  systems  of  die  quarter*] 


mass  to  the  unspaing  mass  is  positive  (i.e.,  damper  extending).  If  the  unsprung  mass 
happens  to  be  moving  upwards  faster  than  the  sprung  mass,  the  damper  can  only  apply 
forces  in  the  wrong  direction. 

Thus,  the  damping  force  (Fd)  desired  of  the  ER  damper  is  proportional  to  Vs 
whenever  Vs  and  V(s.u)  ue  of  the  same  polarity,  but  this  force  is  to  be  zero  whenever  Vs 
and  V(s.u)  are  of  opposite  polarity.  In  equation  form; 

Fd  ss  CgkyVs  •  [(Vs  •  V(s.u)  >  0]  (14) 

where  Cgky  is  a  constant  and  the  logical  function  [(Vs  ■  V(s.U)  >  0]  s  i  if  TRUE  and  0  if 
FALSE. 

Equation  14  indicates  the  desired  damping  force,  but  not  the  ER-damper-command 
voltage.  The  previous  discussion  of  the  static  tests  of  the  damper  showed  the  relationsliip 
between  damper  force  and  voltage  to  be: 

IFdl-lkvvPl  (15) 

where  ky  and  |3  are  both  constants.  Combining  equations  14  and  IS  and  solving  for  the 
voltage  leads  to  the  following  expression  for  the  voltage  command  sent  to  the  high- 
voltage  amplifier: 

Vc  =  lk«Vsl»^*[(Vs'V(s.U)>0]  (16) 

where  k  is  a  constant  equal  to  Cgky/ky. 

Equation  16  is  the  skyhook  control  algorithm  for  the  ER  damper  u.sed  in  our  testing. 
From  the  results  of  static  testing,  the  chosen  value  of  p  was  2.5.  By  preliminary 
experimentation,  15x10®  v^*®  (in/sec)‘l  was  chosen  as  the  value  for  k. 

Figure  13  is  a  plot  of  an  actual  voltage  command  time  history  generated  by  this 
skyhook  control  algorithm  during  a  test  run  in  which  the  quarter-HMMWV  was  excited 
with  a  1.3  hz  sinusoidal  road  motion.  The  sinusoidal-like  quality  of  the  signal  is  apparent 
as  is  the  on/off  switching  function,  which  is  based  on  the  relative  velocity. 

The  conventional  shock  absorber  algoriihm  employed  only  the  relative  velocity  V(s- 
U).  This  algorithm  was  simply  a  representation  of  the  convenlional  bilinear  behavior  of 
automotive  shock  absorbers'  with  the  addition  of  the  non-linear  quality  of  the  ER  damper 
force-to-voltage  relationship. 

The  equations  used  to  simulate  a  conventional  shock  were  simply; 

Vc  =  lkc*V(s-u)l’^P  whenV(s.u)<0  (17) 
Vc  =  Ikt  •  V(s.U)l  when  V(s.U)  >*  0 

where  kc  and  kt  are  constants  for  compression  and  extension,  respectively.  The  constant 
values  actually  used  in  testing  were  kc  =  8x10®  v2  5  (in/sec)'*  and  ki  »  16x10®  v^  ® 
(in/sec)-'. 


'  Conventional  shock  absorbers  typically  provide  stronger  resistance  to  extension  than  to  compression 
(basically  because  hitting  curbs  can  generate  much  higher  unsprung  moss  velocities  than  falling  into  holes). 
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Figure  13.  Voltage  command  to  the  ER  damper  using  the  skyhook  algorithm 

Before  leaving  this  discussion,  consider  the  source  of  the  velocities  used  in  equations 
16  and  17.  Vs  was  determined  by  integrating  the  signal  provided  by  the  accelerometer 
mounted  on  the  sprung  mass.  (On  the  test  rig,  we  could  measure  this  velocity  directly, 
but  in  the  real  world,  integration  of  acceleration  is  the  most  reasonable  approach.)  The 
frequency  domain  transfer  function  shown  in  the  block  diagram  of  figure  12  indicates  that 
the  acceleration  signal  is  simultaneously  subjected  to  a  high-pass  filter  to  remove  any  DC 
component  and  prevent  drift  of  the  calculated  velocity.  To  obtain  V(s.u)i  the  signal  from 
a  displacement  transducer  placed  in  parallel  with  the  damper  was  differentiated  using  a 
simple  digital  algorithm. 

Data  Acquisition 

Table  2  identifies  the  17  data  signals  collected  and  stored  continuously  during  testing 
on  the  dynamic  facility.  The  table  corresponds  to  many  of  the  signals  identified  earlier  in 
figure  12.  Signals  0  through  11  in  the  table  are  derived  directly  from  the  instruments 
listed  in  the  table  as  the  Signal  Source.  Signals  12  through  16  were  derived  through 
continuous  on  line  linear  or  nonlinear  transformations  of  the  first  12  signals,  and  were 
stored  in  data  files  along  with  the  primary  data. 

As  shown  in  figure  12,  signals  9,  10,  11,  and  16  were  isolated  from  the  computer 
system  using  custom  designed  high-voltage  isolation  amplifiers.  These  signals  were 
subject  to  high-voltage  spikes  in  the  event  of  arcing  within  the  ER  shock  absorber.  The 
special  isolation  was  necessary  to  protect  the  electronic  circuits. 

The  high-voltage  amplifier  of  figure  12  was  a  Trek  Model  640.  This  amplifier  is 
capable  of  following  a  sinusoidal  command  input  of  frequencies  up  to  10  kilohz.  h  can 
provide  200  watts  of  output  power  at  up  to  ±10  kilovolts. 


Table  2.  Instrumentation  and  recorded  data  signals 


No 

Signal  Source 

Signal 

Symbol 

0 

5  klb  load  cell 

Road  force  1 

Fri 

1 

10  klb  load  cell 

Road  force  2 

Fr2 

2 

5  klb  load  cell 

Road  force  3 

Fr3 

3 

LVDT 

Road  vertical  displacement 

dR 

4 

Servo  accelerometer 

Road  vertical  acceleration 

aR 

5 

LVDT 

Damper  displacement 

d(s-U) 

6 

Servo  accelerometer 

Sprung  mass  acceleration 

as 

7 

D/A  filtered  output 

Road  vertical  displacement  command 

Cr 

8 

String  pot 

Sprung  mass  vertical  displacement 

ds 

9 

H.V.  amplifier  output 

Instantaneous  voltage  output 

Vm 

10 

H.  V.  amplifier  output 

Instantaneous  current  output 

im 

11 

Thermocouple  &  amplifier 

Internal  damper  temperature 

T 

12 

Derived  data 

Tire  force 

Ft 

13 

Derived  data 

Digitally  filtered  damper  displacement 

d(S-U)f 

14 

Derived  data 

Damper  velocity 

V(s-u) 

15 

Derived  data 

Sprung  mass  velocity 

Vs 

16 

Derived  data 

H.V.  amplifier  command  voltage 

All  data  signals  were  collected  at  a  sampling  rate  of  1700  hz  through  losv-pass,  SO  hz, 
analog  filters.  The  delay  between  individual  channel  samples  was  10  |isec  at  the  A/D 
converter.  The  high  acquisition  rate  was  required  to  ensure  enough  data  for  effective 
nitering  of  derived  damper  velocities.  Derived  data  were  computed  and  updated  by  the 
controller  at  its  command  frequency,  100  hz.  These  data  were  subsequently  over  sampled 
(by  a  factor  of  17)  when  the  file  was  saved  to  produce  a  uniform  1700  hz  data  rate  for  all 
channels. 

The  data  files  contain  12-bit/sample  data  and  conform  to  a  standardized  ERD  file 
:  tructure  which  includes  measured  amplifier  gains  and  offsets  for  each  channel,  as  well  as 
channel  names,  units,  and  instrumentation  notes.  Data  was  reduced  using  Split  Plus,  a 
program  designed  to  filter,  subsample,  and  perform  linear  transformations  on  ERD  data. 
Ail  channels  of  all  files  were  digitally  filtered  at  170  hz  and  subsampled  to  340  hz  prior  to 
any  analysis  to  reduce  the  number  of  data  points  (by  a  factor  of  S). 

The  following  linear  transformation  was  performed  to  produce  the  tire-vertical-force 
signal. 

Ft  *  (Fri  +  Fr2  +  Fr3)  -  m*aR  (18) 

where  m  is  the  mass  of  the  simulator’s  road  clement  supported  by  the  load  cells.  Since 
this  mass  lies  between  the  load  cell  and  point  of  tire-road  contact,  the  last  term  of 
equation  18  is  necessary  to  remove  that  portion  of  load  cell  force  which  is  expended  in 
accelerating  the  road,  leaving  the  normal  force  at  the  tire/road  interface. 
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The  digitallj'  filtered  damper  displacement,  d(s.U).  was  computed  using  a  64-point 
moving  window  average. 

d(S-U)f  =  ^last  64  d(s.u)  samples  acquired)  (19) 

The  relative  time  of  the  last  displacement  sample  acquired  was  noted  for  computation 
of  damper  velocity,  V(s.u)*  Velocity  was  computed  in  in/sec  using  the  equation: 

V(S.U)  *  (  d(S-U)fn  ‘  d(S-U)fn.l  )  /  (^0  (20) 

where  n  is  the  sampling  index  and  At  is  the  time  step  between  data  samples. 

The  sprung  mass  velocity,  Vs,  was  computed  from  the  sprung  mass  acceleration,  as, 
using  the  algorithm  that  combines  integrator  and  high-pass>fllter  functions  and  which  was 
discussed  in  section  3.1.  The  algorithms  used  to  calculate  the  command  voltage  for  the 
high-voltage  amplifier  were  discussed  in  the  previous  section. 

Power  spectral  densities  (PSD)  were  computed  using  Spectrum,  a  program  designed 
by  UMTRI  to  convert  time-sampled  data  to  frequency  domain  forms.  All  road  data  were 
transformed  to  PSDs  covering  the  frequency  range  of  0  to  15  hz  in  steps  of  0.050  hz.  The 
exceptions  are  CHV  01,  CHV  06,  and  the  sinusoidal  tests,  which  could  only  be 
transformed  in  steps  of  0.10  hz  due  to  the  limited  length  of  the  data  files. 

Transfer  function  gains  were  also  found  using  Spectrum.  Gains  were  computed  for 
the  road  acceleration  as  an  input  and  sprung  mass  acceleration  or  tire  force  as  an  output. 
All  gain  function  results  were  subsequently  window  filtered  with  a  window  length  of  0.20 
hz  prior  to  plotting. 

Road  Actuator  Control 

The  road  proflles  were  generated  by  the  same  computer  that  provided  the  data 
collection  service.  The  profile  was  output  as  a  waveform  by  a  hardware-controlled  D/A 
converter.  The  profiles  were  stored  as  columnar  Hies,  each  value  in  a  file  representing  a 
road  vertical  displacement  in  inches.  Flies  could  contain  measured  profiles  of  actual  roads 
or  idealized  profiles  such  as  sine  waves,  ramp  functions,  step  waveforms,  etc.  Vertical 
displacement  values  were  available  at  4-in  intervals  longitudinally  (three  data  points  per 
foot).  The  speed  at  which  values  were  output  to  the  road  servo  depended  on  the  speed  of 
desired  “travel”  over  the  road. 

Vertical  displacement  values  were  converted  from  in  to  integer  values  to  drive  a  12- 
bit  D/A  converter  at  a  gain  of  ±2047/±5  in.  (Values  exceeding  5  in  were  clipped.)  A  self¬ 
recycling  clock  with  a  10  psec  update  time  was  loaded  with  the  proper  timing  count  to 
output  values  at  equal  time  intervals.  At  each  update  interval,  the  next  value  in  the  profile 
was  sent  to  the  D/A  converter.  The  analog  output  of  the  converter  was  filtered  using  a  1 
khz,  single-pole,  low-pass  filter  prior  to  being  passed  to  the  servo  amplifier.  The  filtered 
data  was  redigitized  and  collected  as  the  road  vertical  displacement  command,  Cr. 

Several  preliminary  tests  were  conducted  using  sinusoidal  wave  forms  to  drive  the 
road  simulator.  The  majority  of  testing,  however,  used  road  profiles  provided  by 
TACOM  and  modified  by  UMTRI  to  fit  the  limits  of  the  road  simulator. 
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Figure  14,  Road  profiles  provided  by  TACOM 

Figure  14  shows  the  eight  road  profiles  provided  by  TACOM.  We  understand  the 
road  profiles  are  of  proving  grounds  courses  at  Aberdeen,  Fort  Knox,  and  elsewhere.  The 
eight  profiles  provided  are  designated  by  TACOM  as  follows: 

APQ09  CHVOl  FTKN31  LETS 

APG29  CHV06  FTKN  34 

APG37 

UMTRI  modified  these  road  profiles  simply  by  subjecting  them  to  high-pass  filtering 
as  required  to  constiain  the  total  displacement  to  10  in — the  limit  of  the  road  simulator. 
This  essentially  means  we  removed  the  hills  but  left  the  bumps.  The  start  and  end  of  the 
individual  road  profiles  were  also  trimmed  to  allow  continuous  repetitive  running  of  the 
segment  with  a  reasonable  transition  at  the  end-to-start  interface.  One  example  of  the 
final  form  of  a  profile  (derived  from  the  longest  profile  of  figure  14)  is  shown  in  figure 
IS.  All  of  the  road  profiles  used  in  testing  are  presented  in  appendix  D. 

The  road  simulator  was  not  sufficiently  powerful  to  follow  these  rough  roads  at  any 
speed.  The  limitation  was  essentially  the  flow  capability  of  the  pump  and  was  reflected 
in  the  maximum  simulated  forward  speed.  At  higher  speeds,  the  displacement  command 
would  have  had  to  be  attenuated  in  order  to  stay  within  the  capability  of  the  system.  The 
theoretical  relationship  between  simulated  vehicle  speed  and  the  required  displacement 
gain  to  allow  continual  running  is  show  in  figure  16.  In  practice,  somewhat  better 
continuous  performance  was  attainable,  and  substantial  improvement  could  be  obtained 
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Figure  15.  Modified  road  proflie  FTKN  31 


Table  3.  Road  displacement  gains  and  simulated  velocities  used  in  testing 


Road 

APO 

APO 

APG 

CHV 

CHV 

FTKN 

FTKN 

LET 

Profile 

09 

29 

37 

01 

06 

31 

34 

5 

mmmm 

0.7 

0.7 

0.7 

0.7 

0.7 

0.7 

0.7 

0.7 

18 

13 

33 

36 

57 

28 

26 

17 

Dynamic  Tests  of  the  ER  Damper  Alone 

A  short  series  of  dynamic  tests  were  conducted  on  the  ER  damper  to  determine  its 
coulomb  and  viscous  -friction  response  and  to  characterize  the  response  time  of  the 
controlled  force  to  applied  voltage. 

Figure  17  is  a  photograph  of  the  test  arrangement.  The  servo  actuator  used  for  the 
road  simulator  and  the  ER  damper  were  arranged  in  line  on  the  bed  plate.  The  actuator 
was  used  to  drive  the  damper  directly  through  a  load  cell,  which  measured  the  damper  rod 
force.  An  LVDT  was  used  to  measure  the  damper  sti*oke. 

Tests  were  conducted  with  three  voltage  wave  forms:  (1)  Vc  ■  0  v  for  a  complete  run; 
(2)  Vc  ■  500  V  for  a  complete  run;  (3)  Vc  ®  500  v  at  t  ■  0,  stepping  up  to  2,500  volts  at  one 
third  of  the  compression  stroke  and  returning  to  0  volts  at  two  thirds  of  the  compression 
stoke.  For  each  wave  form,  four  tests  were  run  with  the  stroke  speeds  of  2, 4, 6  and  8  in/s, 
respectively.  Data  were  used  from  the  compression  stroke  only  because  the  special 
plastic  fittings  used  to  isolate  the  ER  shock  electrically  from  the  load  cell  were  not  strong 
enough  in  tension  to  resist  the  larger  loads. 

Tests  conducted  with  the  first  wave  form  (Vc  ■  0  v)  were  Intended  to  characterize  the 
viscous  behavior  of  the  shock.  The  intent  was  to  measure  the  steady-state  response  of  the 
ER  shock  at  the  four  velocity  levels.  Test  results  showed,  however,  that  fairly  large 
mechanical  vibrations  existed  throughout  the  stroke.  An  example  of  the  time  history  of 
shock  force  during  one  such  test  is  shown  in  figure  18.  After  some  informal  experiments, 
the  second  test  series  using  a  command  voltage  of  500  volts  was  added  with  the  hope  of 
providing  less  noisy  data.  This  was  not  particularly  successful.  In  the  end,  steady  state 
results  were  derived  from  both  of  these  test  series  simply  by  processing  the  time  histories 
through  the  appropriate  low-pass  filters.  The  processed  data  corresponding  to  the 
example  of  figure  18  are  shown  in  figure  19. 

The  third  test  series  involved  step  changes  in  command  voltage  from  500  volts  to 
2500  volts  and  back  to  0  volts  during  the  stroke.  The  resulting  time  histories  of  shock 
force  were  (1)  low-pass  filtered  to  obtain  the  steady-state  response  at  2500  volts,  and  (2) 
examined  to  characterize  the  response  quality  of  the  ER  shock. 

The  steady-state  results  from  the  three  test  series  are  summarized  in  figure  20.  This 
figure  also  includes  the  appropriate  (i.e.,  at  0,  500,  and  2500  volts)  zero  velocity  data 
shown  previously  in  figure  9.  These  static  value:,  from  the  earlier  tests  are  generally 
below  the  trend  shown  by  the  steady-state  values.  The  static  values  are  the  sum  of  the 
shock  forces  produces  by  coulomb  friction  (at  zero  volts)  and  the  so-called  Bingham 
yield  stress  of  the  ER  fluid. 
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Figure  19.  An  example  of  a  heavily  flitered  time  history  of  ER  shock  force  during  a 

steady  state  test  (4  in/sec,  0  voits) 
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Figure  20.  Summary  of  thesteady^state  test  results 

Among  the  steady>state  data,  the  shock  force  at  2  in/sec  and  0  volts  deviates  from  the 
trend  and  almost  coincides  with  the  shock  force  at  2  in/sec  and  SOO  volts.  If  this  point  is 
removed  from  the  plot,  a  line  fitted  to  the  0-volt  line  will  be  nearly  parallel  to  the  lines 
fitted  to  the  500-  and  2500-volt  data,  respectively. 

The  slope  of  tljese  fitted  lines  is  indicative  of  the  viscous  response  of  the  shock.  The 
viscous  damping  forces  are  substantial,  and  larger  than  we  had  intended  within  the 
original  design.  This  result  is  largely  because  the  base  fluid  finally  employed  was 
transformer  oil.  Transformer  oil  has  a  higher  viscosity  than  kerosene,  which  was 
originally  expected  to  be  used  as  the  base  fluid. 

The  vertical  separation  of  the  lines  in  figure  20  is  related  to  the  ER  behavior  of  the 
shock.  The  shock  forces  at  2,500  volts  are  roughly  two  times  higher  than  the  zero-field 
forces.  The  shock  forces  at  4,500  volts,  if  achievable  without  arcing,  would  be 
approximately  4  times  those  at  2,500  volts  (i.e.,  approximately  4,000  pounds)  according 
to  the  2.5*^  power  law  (see  equation  15). 

The  dynamic  response  of  the  ER  shock  absorber  is  illustrated  in  the  time  histories  of 
figures  21  and  22.  Each  figure  includes  the  time  history  of  control  voltage  and  current 
and  the  resulting  shock  force.  Figure  21  shows  the  portion  of  a  single  test  when 
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Figure  21.  Timing  of  the  ER  shock  control  current  and  voltage,  and  shock  force 
upon  a  step  change  of  voltage  command  from  500  volts  to  2500  volts 
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voltagt  stepped  from  500  to  2500  volts.  Figure  22  is  from  the  same  test  at  the  time  when 
the  con\mand  voltage  was  stepped  back  down  to  0  volts. 

Theie  time  histories  seem  to  show  that  when  the  step  change  in  command  voltage  is 
requested,  the  change  in  current  flow  is  very  rapid,  but  the  change  in  actual  voltage  lags 
slightly  due  to  the  need  to  charge  the  plates  inside  the  ER  shock.  The  force  appears  to 
develop  virtually  in  phase  with  control  voltage.  Thus,  it  appears  that  the  dynamic 
response  of  this  system  is  limited  by  the  response  of  the  high-voltage  amplifier,  and  the 
force  responds  to  the  voltage  without  appreciable  lag. 

Quarter-HMMWV  Road  Simulator  Testing 

The  principal  test  aeries  conducted  on  the  road  simulator  involved  analysis  of  the  ride 
performance  of  the  quarter-HMMWV  over  the  road  profiles  provide  by  TACOM. 
However,  a  few  preliminary  tests  to  determine  some  fundamental  system  properties  were 
done  with  pure  sinusoidal  road  displacements. 

The  first  of  the  sinusoidal  tests  were  conducted  with  no  damper  installed.  These  tests 
were  done  at  essentially  the  largest  amplitude  of  road  acceleration  that  did  not  cause  the 
tire  to  leave  the  road  surface  at  the  resonant  frequencies.  (Relatively  strong  nonlinear 
qualities  make  the  results  discussed  somewhat  sensitive  to  magnitude.)  These  tests 
showed  the  sprung  mass  resonant  frequency  to  be  1.5  to  1.6  hz  and  the  unsprung  mass 
resonance  to  be  at  9  to  10  hz.  Linkage  friction  and  tire  damping  resulted  in  the  sprung- 
mass  mode  being  damped  at  about  0.1  of  critical. 

The  ER  damper  was  then  installed  in  the  quarter-HMMWV  and  sinusoidal  tests  in  the 
frequency  range  of  the  sprung  mass  resonance  were  repeated  with  no  control  voltage 
applied.  The  results  showed  the  resonant  frequency  to  be  in  the  1.3  to  1.4  hz  range  and 
the  equivalent  damping  to  be  about  0.25  of  critical.  It  was  thus  apparent  that  the  coulomb 
friction  of  the  ER  damper  was  appreciable  in  this  system. 

In  the  main  series  of  tests,  the  quarter-HMMWV  was  exercised  over  each  of  the  eight 
road  profiles  with  four  different  damper  conditions.  These  conditions  were  (1)  ER 
damper  with  Vc  ■  0;  ER  damper  with  the  skyhook  control  algorithm;  ER  damper  with  the 
conventional  shock  absorber  control  algorithm;  and  the  standard  HMMWV  shock 
absorber. 

The  results  of  all  of  these  tests  were  analyzed  in  the  frequency  domain.  The  transfer 
function  between  vertical  acceleration  of  the  road  surface  and  vertical  acceleration  of  the 
sprung  mass  was  determined  for  each  test  condition.  Additional  analyses  were  performed 
to  examine  the  transfer  function  between  vertical  acceleration  of  the  road  surface  and  tlre- 
to-road  normal  force. 

Figure  23  presents  an  example  of  the  spectral  analysis  process  undertaken  for  each 
test  run.  These  data  are  from  the  run  over  surface  APG  29  with  the  ER  dampei  and  zero 
control  voltage.  The  top  graph  in  this  flgure  is  a  plot  of  the  power  spectral  density  (PSD) 
of  the  vertical  acceleration  of  the  road  surface,  that  is,  the  spectral  content  of  the  input  to 
the  quarter-HMMWV.  The  middle  graph  is  a  plot  of  the  PSD  of  the  resulting  vertical 
acceleration  of  the  sprung  mass — the  output  of  the  system.  While  the  spectral  content  of 
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the  input  is  relatively  uniform,  the  content  of  the  output  is  concentrated  in  the  vicinity  of 
l.S  hz.  This,  of  course,  is  because  this  is  the  natural  frequency  of  the  sprung-mass  mode 
of  the  system,  and  the  system  resonates  at  this  frequency.  The  bottom  graph  is  a  plot  of 
the  transfer  function  of  the  system.  This  plot  is  more-or-less  (but  not  exactly),  the  square 
root  of  the  PSD  of  the  output  divided  by  the  PSD  of  the  input,  evaluated  at  each 
frequency.  The  plot  is,  in  fact,  the  gain  of  the  system  as  a  function  of  frequency.  This 
plot  shows  again  that  this  system  resonates,  or  amplifies  the  input,  in  the  vicinity  of  the 
natural  frequency  of  the  sprung  mass.  At  this  frequency,  the  gain  is  greater  than  unity. 
However,  the  system  effectively  attenuates  the  input  at  frequencies  above  the  sprung 
mass  resonance.  Note  also,  that  there  is  a  very  mild  rise  in  the  gain  at  the  unsprung  mass 
resonance  —  around  9  hz.  (This  effect  is  smdl  here,  and  might  appear  insignificant,  but 
a  similar  rise  of  one  amount  or  another  is  present  around  9  hz  in  virtually  all  of  the 
frequency  domain  data  plots  produced  in  this  study.) 

An  ideal  gain  function  for  good  ride  quality  would  be  one  that  strongly  attenuates  the 
dynamic  response  of  the  sprung  mass  at  all  frequencies.  Practically,  the  goal  of  improved 
ride  is  to  reduce  the  gain  appreciably  at  the  sprung-mass  resonance  without  substantially 
increasing  gain  at  the  unsprung-mass  resonance. 

Figure  24  shows  the  system  gain  (similar  to  the  bottom  plot  of  figure  23)  measured  on 
APQ  29  with  three  different  damper  conditions.  The  upper  plot  shows  the  gain  with  the 
ER  damper  installed  but  with  no  voltage  command.  The  center  plot  shows  the  system 
gain  with  the  ER  damper  operating  as  a  skyhook  damper.  The  bottom  plot  is  for  the 
system  with  the  standard  HMMWV  shock  absorber  installed.  These  data  clearly  show 
the  ER  damper  in  the  skyhook  control  mode  to  be  more  effective  at  reducing  the 
amplification  of  acceleration  at  the  sprung  mass  resonance  than  the  standard  shock 
absorber.  A  small  increase  in  response  for  the  system  with  the  skyhook  damper  is  noted 
in  the  mid  range  frequencies  from  5  to  10  hz. 

Figures  25  and  26  address  two  other  ancillary  issues,  again  with  example  data  from 
APO  29.  Figure  25  compares  the  system  gain  with  the  standard  shock  absorber  (bottom 
plot)  to  the  gain  when  the  ER  damper  is  operated  with  the  conventional  shock  absorber 
control  algorithm.  It  is  apparent  that  with  the  control  constants  chosen  (kc  ^  8x10^ 
v[in/sec]'0-^  and  kt  *  16x10®  vlin/secJ-O'*) ,  the  two  shocks  are  quite  similar.  The  point  to 
be  made  is  that  the  control  constant  used  in  the  skyhook  algorithm  (ka  15x10® 
v[in/sec]'^‘^  )  is  of  a  similar  magnitude  and  will  generate  forces  of  generally  the  same  or 
smaller^  magnitude  then  the  conventional  shock  absorber.  The  improved  performance  of 
the  skyhook  approach  can  thus  be  seen  to  result  from  the  strategic  advantage  of  the 
semiactive  control  approach. 

Figure  26  is  similar  to  figure  24,  but  the  gains  shown  ore  of  tire-to-road  normal  force 
with  respect  to  road  acceleration.  This  gain  is  of  interest  because  low  levels  of  dynamic 
components  of  the  tirc-to-road  contact  force  are  desirable  for  enhanced  road  holding. 


2  Sky  hook  damper  forces  can  be  expected  to  be  smaller  than  forces  developed  by  a  conventional 
shock  absorber  with  similar  properties  because  the  applied  velocities  are  smaller.  The  sky  hook  damper  is 
subject  only  to  sprung  mass  velocities  The  conventional  damper  is  also  subject  to  the  velocity  of  the 
unsprung  mass  which  moves  faster. 
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Figure  24.  Example  sprung  mass  acceleration  gains  measured  with  three  OifTerent 

shock  absorber  conditions  for  APG  29 
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Figure  25.  Example  sprung  mass  acceleration  gains  measured  with 
simulated  and  actual  conventional  shock  absorber  for  APG  29 


Although  the  magnitudes  of  this  gain  are  quite  different  from  those  of  figure  24,  the 
relative  performance  of  the  three  types  of  shocks  is  similar  to  that  which  prevailed  in  the 
earlier  figure. 

The  example  data  presented  here  all  came  from  test  runs  using  the  APG  29  road 
profile.  Measured  gain  functions  for  all  of  the  road  profiles  are  presented  in  appendix  D. 
If  the  vehicle  systems  were  linear,  the  gain  functions  would  all  appear  the  same 
regardless  of  the  road,  However,  this  system  contains  substantial  non-linear  elements. 
We  have  seen  that  coulomb  friction  in  the  linkages  and  the  ER  damper  is  significant.  All 
of  the  shock  absorbers  have  significant  nonlinear  qualities.  The  conventional  shocks 
(standard  shock  or  the  conventional  algorithm  applied  to  the  ER  damper)  are  nonlinear 
due  to  different  compression  and  extension  rates.  The  logical  switching  function  of  the 
skyhook  algorithm  is  very  nonlinear.  Finally,  the  large  displacement  road  profiles  cause 
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the  tire  to  operate  in  its  nonlinear  range  (including  tire  lift).  Thus,  the  results  from 
different  roads  are  indeed  different,  but  they  virtually  all  maintain  the  same  relative 
performance  relationships  observed  in  the  examples  presented  here.  Thus,  we  conclude 
that  this  testing  program  has  successfully  demonstrated  the  potential  effectiveness  of  ER 
fluid  technology  in  intelligent  suspension  control  applications. 
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7.  Quarter-Car  Model  of  the  System 

For  this  study,  the  semiactive  system  was  simulated  using  the  quarter  model  as  shown 
in  figure  7.4.  The  complete  simulation  is  provided  as  a  computer  program  listing  in 
appendix  C.  A  review  of  the  simulation  and  brief  discussion  of  calculation  results  is 
presented  in  this  section. 

The  momentum  equation  for  the  sprung  mass  is  simplified  ns 

(21) 

where  is  the  sprung  mass,  the  spring  constant,  the  damping  force  (positive 
during  compression),  t  the  time,  the  sprung  mass  displacement,  Xy  the  unsprung 
mass  displacement. 

The  momentum  equation  for  the  unsprung  mass  is 

Mu-^^~Fd~‘^s(xu-xs)  +  K'j'(xQ-Xij)  (22) 

where  My  is  the  unsprung  mass,  Kf  the  tire  spring  constant,  and  Xq  the  ground 
displacement. 

In  a  shear-mode  damper,  the  influence  of  the  inertia  force  is  relatively  small,  and  the 
damping  force  in  the  strut  F^  can  be  approximated  with  a  Newtonian  damping  force  FdN 

and  an  ER  damping  force  F^£^,  corresponding  to  the  effects  of  the  plastic  viscosity  (r|) 

and  the  yield  shear  stress  (Xy),  itspectively,  in  the  Bingham  model  (Lou  et  al^  1992). 
(See  equations  4  and  S  in  chapter  5.)  The  Newtonian  and  ER  damping  forces  here  are 
completely  uncoupled  and  additive.  However,  a  complete  solution  of  the  problem 
demands  a  full  solution  of  the  fluid  flow  problem  (Lou  et  ai,  1993). 

From  equations  14,  IS,  and  16  (chapter  6),  one  can  derive  that  the  damping 
coefficient  for  the  ER  damper  is  equal  to  k^k.  With  k  s  15x10®  vo//^‘®(i«/sec)  = 
590.55x10®  volts^'^ (m/ see)  and  ky  =  4.445x10®  the  ER  damping 

coefficient  is  2,625  N/(m/s).  If  one  uses  a  sprung  mass  of  998  kg  and  a  spring  constant  of 
0. 1633  X 10®  N/m,  the  ER  damping  ratio  is  0. 1028. 

The  viscosity  of  the  ER  fluid  at  zero-field  was  not  directly  measured.  From  the  slope 
of  the  0-volt  line  in  figure  20,  it  was  caiculated  that  the  zero-field  viscosity  rj  =  222  cP  or 
0.222  Ns/m. 

Similarly,  the  viscosity  variation  of  the  ER  fluid  as  a  function  of  the  electric  field  was 
not  directly  calibrated  in  this  study.  The  value  of  j8  in  equation  15  which  was  2.5,  is 
based  on  the  data  in  figure  9.  These  data  show  only  the  static  response  of  the  ER  damper. 
The  steady-state  data  in  figure  20  may  better  represent  the  ER  effect.  The  related 


information  from  figure  20  is  listed  in  table  4.  The  relationship  between  the  yield  stress 
and  the  electric  field  is  almost  linear.  The  disagreement  between  the  static  and  steady- 
state  data  are  thus  significant.  Both  cases  may  contain  errors.  ER  fluids  are  known  to 
have  erratic  behaviors  at  zero  or  low  shear  rates.  The  existence  of  the  static  yield  shear 
stress  is  the  subject  of  scientific  debate.  Only  two  elevated  voltage  values  are  available  in 
the  steady-state  data,  and  it  is  therefore  difficult  to  rely  on  these  data  to  derive  the  ER 
law.  The  approach  in  the  numerical  simulation  is  a  compromised  one.  We  chose  a  P 
value  of  2,  which  falls  between  1.0  (steady-state  case)  and  2,5  (static  case).  We  have 

Ty-95xl0'^v2  (23) 

where  V  is  the  voltage,  and  the  constant  95  x  10**^  is  derived  based  on  the  yield  stress  at 
2,500  volts.  Data  at  500  volts  were  ignored  because  of  the  proximity  between  500  and  0 
volts  and  the  dominance  of  the  ER  behavior  at  2,500  volts. 

_ Table  4.  Information  from  figure  20 _ 


voltage 

ER  damping  force 

ER  yield  stress  Xy 

500 

2500 

591  N  (133  Ibf) 

2,504  N  (563  Ibf) 

140  Pa  (0.020  psi) 

594  Pa  (0.086  psi) 

The  steady-state  data  were  obtained  and  analyzed  after  the  dynamic  tests  of  the  ER 
damper,  where  the  2.5^^  power  law  was  used  in  controlling  the  voltage  applied  to  the 
damper.  In  order  to  compare  the  numerical  results  with  the  dynamic  test  results,  the 
voltage  in  the  simulation  will  be  controlled  in  the  same  manner,  that  is,  according  to  the 
2.5*^  power  law.  However,  the  ER  response  to  the  electric  field  is  based  on  the  2"‘l  power 
law,  which  is  believed  to  be  also  governing  the  ER  response  during  the  dynamic  tests. 

During  nie  dynamic  tests,  the  maximum  cutoff  voltage  was  set  3,(X)0  volt,  which  is 
also  used  in  the  numerical  simulation.  Based  on  the  2^^  power  law  and  the  ER  yield  stress 
value  at  2,500  v,  the  ER  yield  stress  at  3,000  volts  is  about  855  Pa  (0.124  psi),  which  is 
far  from  ideal.  With  the  zero-field  viscosity  n  «  222  cP  or  0.222  Ns/m,  the  Newtonian 
shear  stress  reaches  855  Pa  if  the  shear  rate  is  at  3,851  s'^ 

Although  ER  fluids  respond  to  the  field  variation  rapidly,  the  response  is  not  instant. 
The  voltage  across  a  capacitive  ER  damper  has  a  delay  of  its  own.  These  two  delays  are 
simply  lumped  together  in  a  delay  time  constant  r^gigy  in  the  following  equation, 

^  =  — (95x10"*V^-tJ  (24) 

^delay 

The  ER  damper  is  simulated  in  any  one  of  the  following  three  control  types,  as 
described  in  chapter  5, 

(1)  Passive  damper  (i.e.,  no  ER  control).  The  electric  field  is  set  at  zero,  and  only 
Newtonian  viscosity  is  left  to  damp  the  relative  motion. 


(2)  The  skyhook  damper.  The  electric  field  is  regulated  to  damp  the  absolute  vertical 
velocity  of  the  sprung  mass, 

(3)  Simulated  conventional  damper.  The  electric  field  is  regulated  in  such  a  way  that 
the  damper  behaves  like  a  conventional  bilinear  shock  absorber. 

A  computer  code  is  developed  using  the  mathematical  model  and  control  strategies 
discussed  above,  The  momentum  equations  are  integrated  over  the  time  domain. 

At  time  zero,  both  the  sprung  and  unsprung  masses  are  set  motionless,  and  the  ground 
displacement  is  zero.  A  sinusoidal  vertical  displacement  is  then  imposed  at  the  ground, 
and  the  simulation  continues  with  many  cycles  until  a  stable  periodic  motion  is  achieved 
by  the  sprung  mass.  More  specifically,  the  simulation  is  stable  when  the  normalized 
deviation  between  the  sprung  mass  acceleration  amplitudes  of  two  consecutive  cycles  is 
within  0.01  percent.  For  most  runs  of  the  simulation,  stable  solutions  are  obtained  within 
10  cycles.  Simulation  output  from  the  last  cycle  is  used  for  portraying  the  results. 

Each  sinusoidal  cycle  is  divided  into  a  finite  number  of  time  steps  of  equal  step-size. 
For  runs  presented  in  this  study,  there  are  360  steps  in  each  cycle,  with  the  step-size  being 
one  degree  out  of  a  360-degree  cycle,  At  each  time-step,  the  numerical  integration  is 
implemented  with  an  explicit  finite-difference  scheme.  The  accelerations  of  the  sprung 
and  unsprung  masses  are  Iterated  with  a  relaxation  factor,  which  is  equal  to  0.2  in  the  runs 
presented  in  this  report.  The  iteration  stops  when  the  normalized  convergence  of  the 
sprung  mass  acceleration  is  within  0.1  percent. 

The  quarter-car  model  was  simulated  with  each  of  the  three  control  types  (i.e.,  no  ER 
control,  skyhook,  and  simulated  conventional).  For  each  control  type,  a  spectrum  of  25 
frequencies  are  tested,  at  the  values  of  0.1,  0.2,  O.S,  0.8,  1.0,  1.5,  1.8,  1.9,  2.0,  2.5,  3.0, 
4.0,  6.0, 7.0, 8.0, 8.5, 9,0, 9.5, 10, 14, 16, 18, 20,  and  25  hz.  Other  major  parameters  used 
in  the  simulation  are  listed  in  table  5.  There  are  40  hot  electrodes,  i.e.,  40  pairs  of  hot  and 
ground  electrodes  facing  each  other.  The  gap  between  each  pair  of  electrodes  is  1  mm. 
The  outer  radius  of  the  shear  surface  is  equal  to  the  outer  radius  of  the  hot  electrodes,  i.e., 
4.83  cm.  The  inner  radius  of  the  shear  surface  is  equal  to  the  inner  radius  of  the  ground 
electrodes,  that  is,  2.95  cm.  kg/m^.  The  ER  fluid  mass  density  is  chosen  to  be  1,000 
kg/m^.  The  Newtonian  viscosity  is  equal  to  222  cP,  which  is  derived  from  the  steady- 
state  data  as  discussed  earlier.  The  sprung  and  unsprung  masses  and  spring  stiffnesses  are 
effective  values,  i.e.,  their  equivalent  values  for  the  quarter  car  model.  The  tire  spring 

stiffness  is  0.4638  xlO^N/m,  and  damping  in  the  tire  is  neglected.  The  skyhook 
command  constant  is  5.9055x10®  volt^'Ss/m.  For  the  simulated  conventional  damper, 
‘he  command  constants  for  compression  and  extension  are  3.1500x10®  and 
6.3000  X 10®  volt^'^s/m,  respectively.  The  time  delay  constant  is  0.2  ms,  which  is  more  or 
less  arbitrary  but  not  critical  for  the  frequency  range  of  the  simulation.  The  peak  voltage, 
which  the  control  voltage  is  not  allowed  to  exceed,  is  3,000  volts.  The  ball  screw  has  a 
lead  angle  of  17.66  degrees  and  a  pitch  diameter  of  2.54  cm.  Its  frictional  coefficient  is 
0,03,  which,  in  combination  with  other  parameters,  results  in  an  efficiency  of  90  percent 
given  by  the  manufacturer. 


57 


Table  5.  Parameters  used  in  the  simulation 


Parameters 

Values 

Number  of  hot  electrodes  n 

40 

Oap  between  electrodes  h 

1  mm 

Shear  surface  outer  radius  r^ 

4.83  cm 

Shear  surface  inner  radius  r^ 

2.95  cm 

ER  fluid  density  P 

1,000  kg/m^ 

Newtonian  viscosity  T| 

0.222  Ns/m2  (222  cP) 

Effective  sprung  mass  Ms 

998  kg 

Effective  unsprung  mass  My 

181kg 

Effective  spring  stiffness 

0.1633x10®  N/m 

Tire  spring  stiffness  Kf 

0.4638x10®  N/m 

Skyhook  command  constant  k 

5.9055x10®  Volt^'h/m 

Conventional  compression  const  kg 

3.1500x10®  Volt^'^s/m 

Conventional  extension  constant  k. 

6.3000x10®  Volt^-^s/m 

Delay  time  constant 

0.0002  second 

Peak  voltage 

3,000  v 

Ground  wave  amplitude 

1  cm 

Ball  screw  lead  angle  X 

17.66  degrees 

Ball  screw  pitch  diameter  Dp 

2.54  cm 

Ball  screw  frictional  coefficient  / 

0.033,  along  with  other  parameters, 
to  give  an  efficiency  of  90  percent. 

Time  histories  of  the  sprung  mass  velocity  the  sprung-unsprung  mass  veiocity 
difference  V^s-u)^  Newtonian  damping  force  F^f^,  and  ER  damping  force  at  an 
excitation  frequency  of  1 .5  hz  and  under  the  skyhook  control  strategy  are  presented  in 
figure  27.  For  the  ease  of  presentation,  Vs  and  V^s-U)  normalized  by  the  amplitude  of 
the  ground  velocity,  and  ^dN  and  are  scaled  down  by  a  factor  of  1000.  The 
Newtonian  damping  force  responds  linearly  (with  an  opposite  polarity)  to  the 
relative  velocity  V^s-u)-  The  ER  damping  force  is  linearly  proportional  to  the 
negative  sprung  mass  velocity  (-  when  Vs  and  V^s-u)  of  the  same  polarity 
(approximately  between  0  and  120  degrees  and  between  180  and  300  degrees  in  figure 
27)  and  equal  to  zero  when  V5  and  V^^s-U)  of  the  opposite  polarities  (approximately 
between  120  and  180  degrees  and  between  300  and  360  degrees  in  figure  27).  For  this 


— o —  normalized  Vg  - /1 000,  N 

»  normalized  Vu)  -o- FdER /lOOO.N 


Phase  Angle,  deg. 

Figure  27.  Time  histories  of  the  sprung  mass  velocity  V5,  the  sprung-unsprung  mass 
velocity  difference  V(^s-U)>  Newtonian  damping  force  and  ER  damping  force 
at  an  excitation  ft-equency  of  1..*  hz  and  under  the  skyhook  control  strategy 

particular  case,  the  peak  ER  damping  force  is  about  twice  the  peak  Newtonian  damping 
force. 

Time  histories  of  the  sprung  mass  displacement,  velocity,  and  acceleration  and  the 
unsprung  mass  displacement,  velocity,  and  acceleration  for  the  same  run  are  plotted  in 
figure  28.  The  velocities  are  again  normalized  by  the  amplitude  of  the  ground  velocity. 
The  displacements  and  accelerations  are  normalized  by  the  amplitudes  of  the  ground 
displacement  and  acceleration,  respectively.  The  variables  of  the  unsprung  mass  have 
generally  more  high-frequency  components  than  those  of  the  sprung  mass  because  of  its 
small  mass.  The  unsprung  mass  is  more  or  less  “pushed  around”  by  the  sprung  mass.  The 
gain  of  the  sprung  mass,  that  is,  the  sprung  mass  acceleration  normalized  by  the  ground 
acceleration  amplitude,  is  1.393. 

The  gain  of  the  sprung  mass  under  (1)  no-ER  control,  (2)  skyhook  control,  and  (3) 
simulated  conventional  control  arc  presented  in  figure  29.  The  skyhook  control  is 
generally  better,  in  terms  of  reducing  the  gain,  than  the  simulated  convention.il  control, 
which  is  in  turn  better  than  no  ER  control.  In  the  vicinity  of  2  hz,  the  gains  of  no-ER 
control,  skyhook  control,  and  simulated  conventional  control  are  2.7,  1.5  and  2.1, 
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Figure  28.  Time  histories  of  the  sprung  mass  displacement ,  velocity,  and 
acceleration  and  the  unsprung  mass  displacement,  velocity,  and  acceleration 
at  an  excitation  frequency  of  1.5  hz  and  under  the  skyhook  control  strategy 

respectively.  The  corresponding  gains  in  the  experiments  (figures  24  and  25)  are  2.7,  2.0, 
and  2.5,  respectively.  The  numerical  simulation  results  generally  agree  with  the 
experimental  results.  In  particular,  the  gain  values  do  not  show  notable  increases  in  the 
vicinity  of  the  wheel-hop  frequency,  9  hz. 

The  frequency  responses  of  the  unsprung  mass  gain  under  (1)  no  ER  control,  (2) 
skyhook  control,  and  (3)  simulated  conventional  control  are  plotted  in  figure  30.  For  the 
unsprung  mass,  the  gain  at  the  vehicle  resonance  frequency  is  higher  under  the  skyhook 
control  than  under  either  no  ER  control  or  simulated  conventional  shock.  The  unsprung 
mass  is  again  “pushed  around”  by  the  sprung  mass,  whose  damping  is  the  dominant 
priority  of  the  skyhook  control  scheme. 
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Figure  29.  The  fk'equency  responses  of  the  sprung  mass  gain  under  (1)  no  ER 
control,  (2)  skyhook  control,  and  (3)  simulated  conventional  control 


I 


Figure  30.  The  frequency  responses  of  the  unsprung  mass  gain  under  (1)  iiu  ER 
control,  (2)  skyhook  control,  and  (3)  simulated  conventional  control 
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8.  Conclusions 


The  working  prototype  of  an  electrorheological  damper  has  been  developed  and 
demonstrated.  This  device  has  shown  a  level  of  suspension  control  that  is  superior  to  that 
of  conventional  shock-absorber  damping.  This  demonstration  is  significant  primarily 
because  it  illustrates  the  system-level  capability  of  £R  technology.  That  is,  it  is  more  or 
less  straightforward  at  the  present  time  to  configure  an  BR  device  to  control  force  or 
moment  over  a  broad  bandwidth  of  responses  in  certain  classes  of  automotive  devices. 
On  the  other  hand,  the  project  has  also  documented  that  performance  limitations  continue 
to  constrain  the  range  of  achievement  with  ER  devices. 

As  the  project  concludes  in  early  1994,  it  is  also  apparent  that  parallel  advancements 
in  ER  fluid  technology  ore  offering  to  extend  the  fluid  performance  range  and  to  make 
possible  entirely  different  design  approaches  for  configuring  an  automotive  suspension 
damper.  For  example,  an  American  chemical  supplier  is  advertising  ER  fluids  with  1 .0 
pounds  per  square  inch  shear  strength  capability.  The  fluids  are  comprised  of  cross- 
linked  polyurethane  particles  in  siiicone  oil.  If  such  a  fluid  was  found  to  exhibit  other 
suitable  attributes,  it  would  allow  the  design  of  a  purely  rotary  damper  that  would  avoid 
the  relatively  high  complexity  and  weight  of  the  ER  damper  that  was  prototyped  in  this 
study.  A  rotary  damper  of  this  type  could  be  fitted  along  the  centerline  of  the  inboard 
pivot  of  either  the  lower  or  upper  control  arm  of  the  HMMWV  suspension. 

The  test  methodology  shown  in  this  study  is  also  instructive  of  practices  that  are 
suited  to  the  measurement  of  the  performance  of  ER  suspension  controllers.  In  particular, 
the  study  employed  certain  novel  optical  isolators  which  protect  equipment  from  the 
high-frequency  voltage  spikes  which  may  attend  inadvertent  arcing  within  the  ER  device. 
Believing  that  arcing  will  be  inevitable  in  any  developmental  setup  with  ER  fluids 
(operating  at  kilovolt  levels  of  excitation)  such  protection  is  axiomatic.  Also,  the  ability 
to  produce  displacements  of  a  simulated  road  surface  which  follow  either  analytic 
waveforms  or  actual  road  profiles  is  useful  for  combining  clinical  studies  of  damper 
behavior  with  realistic  assessments  of  overall  performance.  A  quarter-car  mockup  such 
as  employed  here  seems  entirely  suited  to  ail  initial  examinations  of  ER  suspension 
dampers. 

Finally,  while  this  study  borrowed  from  the  extensive  literature  on  semiactive 
suspension  control,  it  did  serve  to  illustrate  the  merits  of  the  broad  bandwidth  capability 
of  ER  technology  for  the  implementation  of  one  of  the  more  attractive  semiactive 
suspension  algorithms,  skyhook  damping.  The  very  high  speed  response  of  ER  fluids 
does,  indeed,  offer  a  powerful  mechanism  for  implementing  the  rapid  switching  within 
individual  cycles  which  skyhook  damping  requires. 
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APPENDIX  A: 

An  Analysis  of  Electrorheological  Dampers 

This  appendix  is  a  is  l  reprint  of  a  paper  entitled  “A  preliminary  parametric  study  of 
electrorheological  dampers”  presented  by  Lou,  Ervin  and  Filisko  at  the  1993  summer 
meeting  of  the  American  Society  of  Mechanical  Engineers  in  V)' .tshington,  D.C.  It 
presents  the  analysis  which  shows  the  superiority  of  the  shear  mode  of  ER  fluid  operation 
in  the  automotive  damper  application. 
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ABSTRACT 

In  approaching  the  design  of  an  eleetrorheology*bued,  semi*active  suspension,  the 
eleetrorteologicnl  component  (BR  danqier)  can  be  built  u  either  a  tlow*mode,  shear-mode,  or 
mixed-mode  type  of  daaqper.  The  source  of  dashing  force  in  the  flow-mode  is  exclusively  bom 
flow-induced  mssure  dropacrosi  a  valve,  while  that  in  the  shear-mode  is  purely  from  the  shear 
stress  on  a  lUoing  surface,  ^e  dynamics  of  the  fluid  flow  is  included  in  the  derivation  of  the  zero- 
field  damping  forces.  The  control  effeedveoeas  is  found  to  be  strongly  related  to  the  dynamic 
constant  (which  is  proportional  to  the  square  root  of  the  vibration  frequency)  and,  for  shear-  and 
flow-mode  diunpen,  the  ratio  of  the  piston  area  to  the  ctoss-section  of  the  ER  control  gap.  To 
achieve  the  same  poformance,  a  flow-mode  ER  damper  it  not  as  compact  and  efficient  as  a  shear- 
mode  ER  damper.  With  the  same  ER  damping  force,  a  mixed-mode  damper  la  more  compact  than  a 
shear-mode  diunper.  However,  the  mixed-mMe  damper  does  not  have  as  a  low  zero-field  damping 
force  as  the  shear-mode  damper. 


NOMENCLATURE 
A«  charged  area. 

A^,  Aon*  snd  A«  the  flow-,  muted-,  and  shear-mode  charged  areas. 

Aat  erots-section  area  of  an  ER  valve. 

Ap  cross-section  area  of  a  piston. 

Ai.  A?,  As,  At.  Aiun,  As^nt  Atom,  and  Ai«n  sorae  inteimediate  variables. 
C  the  electric  capacitance  of  a  damper. 

Q  the  electric  current  density  of  an  ER  fluid. 

Ch!  an  inteimediate  variable. 

Dp;  ID  inteimediate  variable. 

Ep  peak  electric  field  iticngth. 

Fem  the  total  zero-field  damping  force  in  a  mixed-coode  damper. 

Fexi  the  ER  damping  force  in  a  flow-mode  damper. 

Ferid  the  ER  damping  force  in  a  mixed-mode  dimpa. 

Fbri  the  ER  damping  force  in  a  ahear-mode  damper. 

Fu  the  zero-field  (upper-plate)  damping  force  in  a  shear-  mode  damper, 
for  the  zero-field  factor  of  a  flow-niDde  damper, 
fan;  the  zero-field  factor  of  a  mixed-mode  dasq>er. 
fwt  the  zero-field  factor  of  a  shear-mode  damper. 
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Si  the  dynaxsic  factor  of  the  lower-plate  stress  of  a  shear-mode  damper. 

gm  die  dynanie  factor  of  Apa 

So  the  dynamic  fhctor  of  a  flow-mode  datx^ 

gom  the  dynamic  factor  of  a  mixed-mode  damper 

gu  the  dynamic  factor  of  (die  upper-plate  rreas  of)  a  shear-mode  damper. 

gun  the  dynamic  factor  of  the  upjw-pw  stress  df  a  mixed-mode  damper. 

U  the  dynamic  factor  of  Apw. 

n  the  gap  size  of  an  ER  valve  or  gap, 
h*  the  dynamic  constant. 

L  the  length  of  an  BR  valve  or  gap. 

P«p  the  peak  capacitive  electric  pcnm. 

Prp  the  peak  restive  electric  power, 

t  the  time. 

Va  charged  volume. 

Vc(i  V«n>  ind  Vti  the  flow-,  mixed-,  and  shear-mode  charged  volumes. 

V  the  piston  velocity. 

Vp  the  amplitude  of  V. 

w  the  widdiofanERvalveorgiq}. 

Xp  the  amplitude  of  the  piston  displaoemeat. 

Apm  the  ER  pressure  drop. 

Apm  the  pressure  drop  component  induced  by  fluid  inertia. 

Ape  the  zero-field  preuure  drop  in  a  flow-mode  damper. 

Api  the  zero-field  pressure  drop  for  an  equivalent  steady  flow  of  the  same  flow  amplimde. 

Apw  the  preuure  drop  component  induced  by  the  wall  shear  streu. 
e  the  dielectric  constant  of  an  ER  fluid. 

^  the  phase  advance  of  the  lower-plate  streu  of  a  shesMnode  damper, 
de  the  phase  advance  of  a  flow-mode  damper. 

0om  the  phase  advance  of  a  mixed-mode  dan^. 

du  the  phue  advance  of  (the  upper-plue  streu  of)  a  shear-mode  damper. 

^  the  phaw  advance  of  the  uppe^plate  streu  of  a  mixed-mode  damper. 

Aw  the  phau  advance  of  Apw. 

T)  the  zero-field  or  pludc  viscosity, 
p  the  density  of  an  ER  fluid. 

X|  the  lower-plate  shear  suest  of  a  shear-mode  daxrqier. 

equivalent  to  the  zero-field  shear  streu  between  two  parallel  plates,  with  one  of  them  fixed 
and  the  other  plate  moving  at  a  constant  velocity  Vp. 
tu  the  upper-plate  shear  stress  of  a  shear-mode  dunpa. 

Xm  the  upper-plate  shear  streu  of  a  mixed-mode  damper. 

Xy  the  yield  suets. 

the  dimensionless  yield  stress. 

(0  the  angular  velocity  of  the  input  vibiadon. 


INTRODUCTION 

An  electrorheologlcal  (ER)  damper  is  a  device  that  is  able  to  generate  a  controllable  damping 
force  by  employing  an  ER  fluid.  An  ER  fluid  can  change  iu  rheological  properties,  most  notably  iu 
effective  viscosity,  within  one  millisecond  when  exwied  to  an  electric  field.  There  have  been 
numerous  designs  of  and  studies  on  ER  dampers  (Brooks,  19g9;  Bullough  and  Foxon,  1978; 
Duclos,  1988:  Hansock  ei  ai,  1991;  Morishiu  and  Miuui,  1992;  Petek.  1992;  Shul'man  et  ol„ 
1987;  Suuiway  er  al.,  1989;  Sianway  et  al.,  1987:  Stevens  er  ai.,  1984).  The  wide  range  and  the 
high  bandwidm  of  the  effective  viscosity  vuiadon  of  ER  fluids  are  two  of  the  primary  factors  that 
have  caused  intense  research  and  development  activities  in  ER  dampen. 

Despite  the  appearance  of  various  new  theological  models  or  ER  fluids  (Oamota  and  Fillsko, 
1991;  Jordan  et  ai.,  1992;  Rajagopal  and  Wineman,  1992),  the  Bingham  model  remains  attractive 


owing  to  its  simplici^  and  effectiveness  (without  significant  loss  in  accuracy)  in  describing 
macroscopic  engineering  problems.  Viscoeiudc  fluid  model  has  also  been  used  in  modeling  ER 
dampen  (Shul'min  ttal„  1987). 

Based  on  the  Bin^am  model,  Stevens  et  al.  (1984)  and  Stanway  et  al,  (1987)  assigned  two 
terms  to  the  das^g  ttiroe  of  an  ER  danger,  with  one  of  them  (a  linear  viscous  damping)  associated 
with  the  plaide  viscosity  and  the  other  (a  Coulomb  friction)  with  the  Bingham  or  yield  stress.  This 
Coulomb  friction  approach  did  not  quite  fit  experimenul  dau,  and  a  power  law  model  wu 
developed  later  (Stanway  et  al,,  1989).  u  the  power  law  model,  the  anqilitude  of  the  damping  force 
it  proportional  to  the  absolute  value  of  the  damper  velocity  to  the  n*^  power.  The  resulu  were  still 
not  satisfaotory  even  with  the  propcttional  coefficient  and  etmonent  being  adjustable. 

The  ER  damper  studies  to  far  have  genensUy  ignored  the  iniluenee  of  the  fluid  dynamics,  i.e., 
the  unsteady  nature  of  the  fluid  velocity,  shear  rates,  and  thus  the  damping  force,  which  are  all 
induced  by  the  inertia  of  the  fluid.  Lou  «  al,  (1992)  have  shown  that  the  fluid  inerda  can  limit  the 
bandwidth  of  an  ER  valve  to  a  level  far  below  thu  the  bandwidth  of  the  rheological  variation  of  an 
ER  fluid.  Similar  influence  exisu  In  ER  dampers,  many  of  which  contain  ER  valves. 

Eleetrorfaeologioal  dainpers  can  be  claatified  as:  flow*mode  (Fig.  la),  mixed-mode  (Fig.  lb), 
and  shear-mode  (Fig.  le).  The  eonfigurationa  illustrated  in  Fig.  1  axe  essentially  generic  concepts 
that  contain  siamlinM  features  to  facilitate  the  design  analyala.  In  pnedoal  hardware,  for  example, 
the  cylinder  la  likelv  to  be  single,  rither  than  double-ended,  and  the  sR  control  valve  or  area  is  likely 
to  be  comprised  of  multiple  parallel  or  concentric  plates.  Among  the  three  conceptual  modes,  the 
flow-mode  damper  is  most  similar  to  the  traditional  shock  absorber  except  that  it  replaces  the 
conventional  orifice  with  an  ER  control  gap  or  valve.  The  source  of  the  damping  force  in  this  mode 
is  exclusively  from  flow-induced  pressure  drop  across  the  piston.  The  ER  valve  can  be  placed 
outside  the  cylinder  (u  in  Hg.  la)  or  within  the  piston  eomponent  For  the  sake  of  compact 
packaging,  the  external  mptoaoh  could  be  ImplementM  by  fcsming  a  flow  chamber  along  the  ouuide 
of  the  cylinder  wall  (Petek,  1992). 

In  the  mixed-mode  dainper,  the  ER  control  is  realised  in  the  gap  between  the  piston  sidewall  and 
the  cylinder  liner.  The  fluid  shear  streu  at  the  shear  surface  also  contributes  to  the  damping  force, 
while  the  flow-lnduosd  pressure  drop  is  still  present  u  a  force  eomponent  on  the  piston. 

If  a  large  by-pus  port  la  cut  through  the  piston  of  a  mixed-mode  damper,  the  device  becomes  a 
simple  shear-mode  element  Since  the  fluid  can  flow  freely  through  the  by-pus,  the  pressure  drop 
across  the  piston  becomes  negligible  and  only  shear  atreaies  developed  along  the  wall  contribute  to 
the  damping  force. 

The  pumse  of  this  smdy  is  to  systematically  analyze  the  dynamic  performances  of  all  three 
modes  of  ER  dampers.  Special  anenuon  is  paid  to  (a)  the  relative  magnltudu  of  the  zero-field  and 
ER  components  ot  the  damping  force  under  dynamic  situation  and  (b)  and  their  implication  in  the 
sizing  and  selecdon  of  the  darnping  mode.  For  simplicity,  ER  fluids  are  usumed  to  be  Bingham 
plastic.  In  cues  where  ER  fluids  are  not  of  Bingham  type,  the  usumpdon  here  will  not  result  in 
significant  etrors  for  p*acticid  dampers  that  far  out-stroke  the  elude  strain  range  and/or  the  low  shear 
rate  range. 


ANALYSIS 
Flow.Mode  Dampers 

If  the  vibradon  of  the  piston  (Fig.  la)  can  be  described  u 

V  ■  VpSin  (Ot  ■>  (oxpsiii  (Ot  ( 1 ) 

where  v  is  the  piston  velocity,  Vp  the  amplitude  of  v,  (o  the  angular  velocity,  t  the  time,  and  Xp  the 
amplitude  of  the  piston  displacement.  From  continuity,  the  amplitude  of  the  (space-wise)  mean 
velocity  of  the  ER  fluid  in  the  ER  valve  is  (Ar/ABx)Xp^  where  Ap  is  the  cross-secdon  area  of  the 
piston,  and  Abr  is  the  cross-secdon  area  of  the  ER  valve.  The  zero-field  (i.e.,  with  electrical  field 
stiengdt  equal  to  zero)  pressure  drop  for  an  equivalent  steady  flow  of  the  same  flow  amplitude  is 


(2) 


where  L  and  h  are  the  length  and  gap  size  of  the  ER  valve,  respecdvely .  The  pbiuneter  t)  is  the  zero- 
field  or  pludc  viscosity,  which  is  usumed  to  be  Newtonian  although  certain  non -Newtonian 
behavior  hu  been  reported  (Lou  et  al,  1990).  This  pressure  drop,  frequently  used  in  the  literature, 
does  not  account  for  the  effect  of  the  dynamics  of  a  Musoldal  flow. 
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Rg.  1.  The  generic  concepa  of  the  three  model  of  electrorheological  dampen. 
A  dynamic  preuuie  drop  hu  a  component  induced  by  fluid  inertia: 

Apm  -  gmAp,lin(0)t  +  J) 


with 


(3) 


(4) 

(5) 


where  p  ii  the  denilty  of  the  ER  fluid.  The  dimeniionleM  piremeter  h*  it,  if  ignoring  the  conitant 
HI,  the  rido  of  the  gap  liie  h  to  the  parameter  group  VTi/pw.  The  lattw  hat  the  dlmeniioni  of 
length.  It  it  related  to  the  location  of  an  ovenhoot  in  the  mean  velocity  of  an  omilla^ 
llow^mode  damper  or  to  the  penetration  depth  of  an  oiclllaioiy  ihear  flow  in  a  ihw-mode 
One  can  thiawnSder  h*  u  a  dimentiooleii  gap  tim.  The  pamneter  h*  and  ia  vai^oni «  caUei 
in  the  literature,  the  Womertley  number  or  dynamic  coniant  or  dlmeniionleii  frequency  becauie  of 

the  exiiance  of  the  frequency  m.  They  are  alio  related  to  the  ratio  of  the  toe 

(which  dependi  on  the  gap  il»  and  the  fluid  Idnematic  viacoilty)  to  the  period  of  the  forcing 

funrton  Ojou  arid  Yang,"  W3).  In  thii  itudy.  h*  U  to  be  called  the  dynamic  constant  instead  of  the 
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dimensionless  gap  size  because  (a)  the  parameter  appears  only  in  a  dynamic  situation,  (b)  the 

practical  gap  size  is  generally  in  a  narrow  range  around  1  mm.  and  (c)  the  frequency  to  changes 
greatly  in  most  ER  fluid  ^tpU^ons. 

The  function  ga.  called  the  dynamic  factor  of  Apn.  is  0  at  h*  ■  0  and  rises  proportionally  with 
h*^  (or  with  to  if  other  parameten  in  h*  is  fixed).  Physically,  the  ratio  of  the  amplitude  of  the  inertia 
component  of  the  dyniuolc  pressure  to  the  equivalent  steady  pressure  drop  is  proportional  to  the 
fieoMnev. 

The  dynamie  pressure  drop  has  another  component  that  is  Induced  by  the  wall  shear  stress; 


with 


gw' 

dw' 

Ch- 

Dh' 


-  gwApiSln(mt  4>  ^) 

(6) 

,2,,.l(2Q,.Ci|-Di|)V4tliil‘" 

3  (2-Ch)*  +  Dg 

(7) 

|+tan-»[21V(2CE-C5.DS)] 

(8) 

_ _ sin  h’eos  h*  4  ainh  h*eosh  h* 

h*[(cosh  h*coi  h*)*  +  (ilnh  h'sln  h*)*] 

(9) 

_ _ sin  h*eoa  h*  -  alnh  h*flQih  h* 

h*[(coih  h*coi  h*)*  +  (linh  h*sin  h*)*l 

(10) 

The  function  gw,  called  the  dynamic  factor  of  Apw.  itaits  from  1.0  at  h*  «  0,  increues  with  h* 
slightly  until  about  h*  ■  1,  and.  then,  rises  almost  prcreoitionally  with  h'*  (or  with  Vtt  If  other 
parameters  in  h*  is  fixed).  Physleally,  the  amplituae  of  the  dynamie  wall  shear  streit  is 
approximately  equal  to  that  of  the  equivalent  steady  shear  atreu  at  low  frequeneiei  and  grows  larger 
athighfpequendei.  -n  . 

The  sexthfield  dynamic  preuure  drop,  the  summation  of  ita  two  components  Apw  and  Apmi  is: 


dpe  ■  geApilin((Dt  -*■  4o) 


with 


3t(2.Ch)*  +  Dii] 


1/3 


do-f  +  ttn-»[IV(2*Ch)] 


(11) 

(12) 

(13) 


where  g«,  and  0e  are  the  dynamic  factor  and  Uie  phue  advance  of  the  flow-mode  damper, 
respectively.  The  primary  value  of  go  is  from  gw  at  low  frequencies  and  from  gm  at  high  frequencies. 
The  factor  g«  remains  close  to  one  at  low  frequencies  and  rises  proportionallv  with  h*^  (Fig.  2). 
Physically,  the  zero-field  dynamic  pressure  drop  is  mainly  to  overcome  the  wall  stress  and  the  fluid 

inertia  at  low  and  high  frequencies,  respectively.  The  phise  sdvance  0e  is  0  at  h*  ■  0  (steady  lUte), 

43”  at  about  h*'  ■  1.  and  90”  as  h*  approaches  infinite. 

The  above  analysis  has  been  for  the  zero-field  situation.  With  the  application  of  an  electric  field, 
an  yield  itreu  is  Induced  in  the  ER  fluid.  For  a  steady  Couette  flow  within  a  narrow  gap,  the  yield 
stress  and  the  aero-field  shear  streu  can  be  superposed.  The  linear  summation  is  no  longer  valid 
under  dynamic  conditions  and  a  precise  solution  can  be  obuined  only  by  solving  the  fluid  dynamic 
equation  fl^u  tt  at.,  1993a).  An^ydcal  solutions  for  a  dynamic  Bingham  fluid  flow  are  not  readily 
available,  while  numerical  solutions  are  generally  dme-consuming  and  not  as  insightful.  For  s 
preliminiry  engineering  parametric  study,  a  simplified,  linear  analytical  approach  is  adopted  here. 

The  ER  pressure  drop  Aper  is  estimated  to  be 


ApER  ■  2^^ 


(14) 
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Fig.  2.  The  dynamic  factor  go  and  the  phaae  advance  of  the  leto-fleld  pieaiure  drop  venui 

die  dynaxtdo  ooutant  b*. 


where  ty  ia  the  yield  atreat.  The  redo  of  the  zexo*field  preaiure  drep  to  the  ER  preaaure  drop  ia 
derived  ai 


■  fef  8in((0t  +  4*o)  (iS) 

Apm 


with 


(16) 

(17) 

(18) 


where  tot  i>  the  zero-field  factor  of  the  flow-mode  damper,  which  ii  a  relative  meuure  of  the  zero- 
field  preaaure  drqp.  For  the  ER  damper  to  be  effective,  the  zero-field  factor  fei  hu  to  be  much  leai 
than  one.  Ihe  new  panmeter  ia  a  dlmensionleu  yield  atreu;  Xo  ia  equivalent  to  the  zero-field 
shear  atreu  tetween  two  puallel  plates,  with  one  of  them  fixed  and  the  other  plate  moving  at  a 
constant  vdodty  Vp. 

The  ER  damping  force  on  the  piston  of  the  flow-mode  damper  is  calculated  u 


FsRf  ■  (19) 

and  it  hu  to  be  able  match  the  designed  peak  damping  force  of  the  damper.  The  compactness  and 
efficiency  of  the  device  are  important  u  well.  The  charged  area  A«  (where  the  charg^  electrodes 
face  each  other)  is  calculated  u 

A«  ■  wL  (20) 


where  w  is  the  width  of  the  ER  valve.  The  charged  volume  Vg  (i.e.,  the  volume  between  the  charged 
electradei)  is  calculated  u 

Ve-hA«-wLh.  (21) 

The  charged  area  Ac  is  related  to  the  electric  c^tacitance  of  the  damper  C  by  the  equation 


(22) 
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where  e  is  the  dielectric  constant  of  the  ER  fluid.  Assuming  a  sinusoidal  electric  field  strength  with  a 
peak  value  of  Ep,  one  obtains  the  following  peak  capacitive  electric  power  Pep  and  peak  resistive 
electric  power  Pip,  respectively, 

Pcp-J«d^Ve  (23) 

and 

Pfp-CdEpVe  (24) 

where  is  the  electric  euirent  density. 

ShBT.MnHa  namnefe 

A  praetioal  shear*mode  damper  hu  a  narrow  gap  whether  it  has  a  concentric  or  parallel 
arrangement  of  itt  electrodes.  The  fluid  flow  in  a  small  gap  can  be  approximated  with  an  one- 
dimensional  flow  between  two  infinite  parallel  plates,  with  one  plate  (denoted  wrbitrarily  as  the 
lower-plate  here)  fixed  and  the  other  plate  (the  upper-plate)  oscillating  in  its  plane  with  a  velocity 

V  ■  VpSin  mt  (25) 

Adapting  the  loludons  firom  the  literature  (Sdiiag  et  al.,  196S)  to  the  current  arrangement  and  a 
Newtonian  fluid  (l.e.  with  zero-field),  the  shear  streaa  at  the  upper-plate  tu  is  derived  as 

(25) 


g«-2VTh'2^^. 

Va^+aJ 

^  ^  +  tan‘*(Aa/Ai)  ■  um''(A4/A3), 

Al  ■  cosh  2h*  cos  2h*, 

A) « sinh  2h'*'  cos  2h*, 


A3«tinh2h*ain2h*, 
Aa  ■  cosh  2h*  tin  2h* 


where  gu  and  00  are  the  dynamic  factor  and  the  phase  advance  of  (the  upper-plate  stress  oO  *he  shear- 
mode  damper,  respectively.  The  dynamic  factor  gu  is  approximately  equal  to  1  when  h*  <  0.5  and 
equal  to  2Vr  h*  when  h*  >  0.5  (Fig.  3).  Therefore,  the  upper-plate  shear  stress  increases  with  the 
dynamic  constant  h*. 


Fig.  3.  The  dynamic  factor  (gu)  and  the  phase  advance  (W  of  the  upper-plate  shear  stress  of  the 

shear-mode  damper. 
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The  lower*plate  shear  stress  X\  is  dehved  as 
ti  -  gite  sin  (cot  -t-  ^ 


(30) 


with 

gi 


2V7h*-sis 


♦i«  J-ian’WAs) 


(31) 

(32) 


where  gj  and  ^  are  the  dynaiaie  factor  and  the  phase  advance  of  the  lower*plate  stress  of  the  shear* 
mode  damper,  respectively.  The  dynamic  factor  gj  ia  qtproximately  equal  to  1  when  h*  <  0.5  and 
approaches  0  at  h{|^  h*  values  (Fig.  4).  Therefore,  the  Iower*plate  shear  stress  decreases  with  the 
dynamic  constant  h*.  During  a  hlgh*fie^eaoy  oadllation,  the  shear  wave  does  not  fully  reaches  the 
lowe^plate,  and  the  tnsiority  of  the  enetgy  ia  oissipated  in  the  fluid  volume  immediately  below  to  the 
upper-plate. 


I 


Fig.  4.  The  dynamic  factor  (gi)  and  the  phase  advance  (^  of  the  lower-plate  shear  stress  of  the 

shw-mode  damper. 


Analytical  solutions  for  dvnamie  flow  of  a  Bingham  fluid  induced  by  an  oscillating  plate  are  not 
available.  Again  u  in  the  analysis  of  the  flow-mode  damper,  a  timplifled.  linear  analyncal  approach 
is  adopted.  In  the  shear-mode  damper,  no  pressure  dm  is  Involved,  and  the  damping  force  is 
directly  related  to  the  shear  force  on  the  upper-plate.  The  ER  damping  force  in  the  shear-mode 
damper  Fgiu  is  esthnaied  as 

Feri  ■  Am^.  (33) 

where  Am  is  the  charged  area  of  the  shear-mode  damper.  The  zero-field  (upper-plate)  damping  force 
is 


Fy  ■  AMtegu  sin(iM  -f  4ii)  (34) 

The  ER  damping  force  Fbri  has  to  be  able  to  match  the  desired  damping  force  of  the  system,  and 
the  zero-field  damping  force  Fy  has  to  be  relatively  small.  The  ratio  of  the  zero-field  damping  force  to 
the  ER  damping  force  is 

■  f;,  sin(tot  •♦■  ^)  (35) 

rgRi 


with 

foi 


ISL 


(36) 
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where  foi  is  the  zero-field  factor  of  the  shear-mode  damper.  Far  the  damper  to  be  effective,  foi  has  to 
be  much  less  than  one. 

Shear-mode  dampen  can  be  further  divided  into  rotary  and  translational  types,  either  of  which 
can  have  in  own  vaiiuons  in  the  way  electrodes  are  arranged. 

The  damper  in  Pis.  Ic  is  only  one  possible  version  of  the  translational  type  of  shear-mode 
dampers,  where  the  elMttodes  are  arranged  concentrically.  Electrodes  can  also  be  parallel  plates. 
Multiple  electrodes  are  generally  used  to  achieve  a  compact  design.  Because  of  the  translational 
motion,  the  cylinder  in  Fig.  Ic  hu  to  be  longer  than  ±e  piston  to  keep  a  constant  charged  area.  This 
results  in  a  significant  amount  of  idling  surface  on  die  cylinder  at  any  moment. 

A  rotary-type  shear-mode  damper  Is  able  to  reduce  die  extra  surface  area  by  having  all  the  surface 
area  engaged  all  the  time.  One  possible  version  it  to  have  a  plurality  of  roudonal  shearing  disks 
along  the  roudon  axis  (L«u  tt  1993b).  An  alternative  is  to  have  multiple  concentric  electrodes 
around  the  roudon  axis.  A  trantladon-to-roudon  device  is  needed  if  the  modon  to  be  damped  is 
transladonaL 


Because  of  the  linear  nature  of  the  flow  of  a  Newtonian  fluid  between  two  parallel  plates,  the 
velocity  between  the  eleetrodes  of  a  mixed*mode  damper  can  be  obtained  by  auperposing  those  of 
the  corresponding  shear-  and  flow-mode  dampen.  The  shear  stteu  on  the  upper-plate  of  a  mixed¬ 
mode  damper  Xw,  it  the  tummadon  of  the  shear  stresses  induced  by  Couette  and  PoiaeuiUe  flows; 


Xim-folwisinfmta-^ 


(37) 


with 


**■ 

0um*tan-‘^ 
Alum 


(38) 

(39) 


Aiim«gu0os9u'^6^gwcot^  and  Ai^n "  Bu  tin  ^  6^^wtin 


(40) 


where  gun  and  Put)  are  the  dynamic  factor  and  the  phase  advance  of  tun,  respecdvely. 
The  total  zero-field  damping  force  Fon  ia  derived  as 


Font  ■  ApgoApiSin  (COt  ■«■  ^)  4  A«gumtoSln  ((«  4  ^ 
and  the  £R  control  force  Fgun,  Is  as  follows 

pBiun  ■  2^Apty  4  AcXy 


(41) 


(42) 


In  either  of  equadons  (41)  and  (42),  the  first  term  on  the  right  hand  tide  is  the  force  from  the 
pressure  drop,  while  the  second  term  represent  the  fence  on  the  theaiing  eurfece. 

The  comblnadon  of  equtdont  (41)  and  (42)  gives 


FeRih 


I  fomSlnfOJt  4  ^ub) 


(43) 


with 


Va?.».4  A? 


(44) 


• -ism  r^ton 
12(Ap/ABii)4l] 


bo 


I  tan 


•  1  Ainm 


Aje 


(45) 

(46) 


Atom  -  12(Ap/AER)*go  COS  0o  +  gwnCOS  bum 
Aiom  -  12(Ap/AER)*Eo  sin  do  4  gumSin  bun 


(47) 

(48) 


Fig.  5.  (a)  The  dynamic  factor  gam  and  (b)  the  phaM  advance  0em  of  a  mixed>inode  damper 
venui  the  dynamic  cooitant  h*  and  the  area  ratio  Aj/Aen* 


In  active  damping  control,  one  likes  to  generate  a  large  ER  damping  foice  and,  at  the  same  time, 
keep  the  zero-field  dampinf  force  as  small  as  possible.  Fm  the  same  talk,  a  shear-mode  damper  and 
a  flow-mode  damp^  therefore  should  have  (a)  the  same  ER  damping  force  and  (b)  the  same  zero- 
field  faaor  value.  For  the  second  requirement,  one  hu 

(49) 

Xoi 

The  substitution  of  equations  (IS)  and  (36)  into  equation  (49)  gives 


An  ■  ^Abr. 
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To  meet  the  first  requirement,  one  can  derive  from  equations  (19),  (20).  (33)  and  (50)  the 
following 

Vc  A«  lu  ^51) 

where  Vm  is  the  charged  volume  of  the  shear-mode  damper,  and  Act  the  charged  area  of  the  How- 
mode  damper.  The  both  dampen  are  assumed  to  have  the  same  gap  size.  The  redo  Vef/Vei  is  a  good 
indicator  of  relative  design  meiltt  of  the  two  dampen  because  the  charged  volume  is  directly  related 
to  the  size  (Fig.  Id)  and  power  requirement  of  a  damper  as  shown  in  equations  (23)  and  (24).  The 
ndo  VrfA'«,  like  go  or  gu,  is  a  function  of  the  dynamic  constant  h*  (see  Fig.  6).  The  ndo  VefA'e$  is 
3  at  h*  ■  0,  decreases  gndually  to  iu  minimum  value  of  1.35  around  h*  -  1,  and  then  increases 
linearly  with  h*.  To  achieve  the  same  pertonnanee,  the  flow-mode  ER  damper  is  therefore  not  as 
compact  and  efficient  as  the  shear-mode  ER  damper. 


Fig.  6.  The  ndo  of  the  flow-mode  ehaxged  volume  (V(<)  to  the  shear-mode  charged  volume 
(Vm)  u  a  funedon  of  the  dynamic  conatant  h*. 

The  slgnlflcanoe  of  the  dynamie  constant  h*  is  evident  in  Fig.  6.  A  praedeal  range  of  h*  it 
etdmated  using  equation  (S)  and  plotted  in  Rg.  7,  where  the  fluid  density  p  and  the  gap  size  h  are 
chosen  to  be  1000  kg/m^  and  1  mm,  respectively.  The  frequency  ranges  from  1  to  1000  Hz,  and  the 
plutic  viscosity  n  ranges  from  10  to  1000  cP  in  Fig.  7.  In  passenger  can,  body  natural  frequencies 
are  1-2  Hz,  the  wheelhop  frequency  is  around  IS  Hz,  and  the  engine  mount  natural  frequency  it  10- 
20  hz  (Morishiu  and  wuui,  1992).  One  thousand  hertzs  is  probably  the  upper  llinit  of  natural 
frequencies  of  any  practical  mechanical  system.  One  hundred  centipolies  is  a  middle  range  value  of 
the  zero-field  vlieoilty  of  ER  fluids.  For  the  chosen  density  ind  gap  size,  the  dynamic  constant  h* 

varies  from  0.05  to  S.O.  Practical  values  of  p  and  h  will  not  devicte  too  much  from  those  chosen 
here,  and  the  range  of  h*  pieiented  in  Fig.  7  is  fairly  r  spreienutive. 


I  u  I  rm  ithT*!' n 


For  the  ER  damping  force  of  a  mixed-mode  damper  Pbrid  and  that  of  a  shear-mode  damper  Fbri 
to  be  equal,  one  obtains 

^  .  6m.  - - 1- -  (22) 

V«  Am  2(Ap(AER)  +1 

where  Van  and  A<m  are  the  charged  volume  and  area  of  the  mixed-mode  damper,  respectively.  The 
gap  size  h  is  assumed  to  be  the  same.  With  a  positive  Ap/AEs,  the  ratio  VemA'c*  is  less  than  one. 
With  the  same  ER  damping  force,  the  mixed-m^  damper  ii  therefore  more  compact  than  the  shear- 
mode  damper. 

The  ratio  of  the  zen>-field  factor  of  the  mixed-mode  damper  fom  to  that  of  the  shear-tnode  damper 
foi  Is  equal  to  gom/gu  because  of  equations  (36)  and  (44).  The  ratio  fon/Toi  is  a  function  of  the 


A-13 


dynamic  constant  h*  and  the  area  ntio  Ap^Aek,  and  is  always  greater  than  one  for  a  positive 
(Fig.  8).  The  mixed-mode  damper  therefore  does  not  have  as  a  low  zero-field  damping  force  u  the 
compinble  shear-mode  damper. 


Fig,  7.  The  dynamlo  constant  h*  u  a  function  of  the  vibration  frequency  (-  u/2n)  and  the  zero-field 
viscosity  (Ti),  with  the  ER  control  gap  size  h  ■  1  mm  and  the  ER  fluid  density  p  «  1000  kg/m^. 


Fig.  8.  The  redo  of  the  zero-field  faaor  of  the  mixed-mode  damper  fem  to  that  of  the  shear-mode 
damper  fgi  as  a  function  of  the  dynamic  constant  h*  and  the  area  redo  Ap/Aot. 


DISCUSSIONS 

Theptimai:y  message  from  the  above  analysis  it;  to  avoid  the  flow-mode  and  to  use  the  shear- 
mode.  There  are  two  physical  reasons  for  the  metuge.  First,  with  a  limited  ER  induced  increue  in 
the  shear  atreu  (i.e.,  Xy),  a  device  ihould  operate  at  as  a  low  shear  rate  as  possible  to  reduce  the 
zero-field  damping  force  and  increue  its  effective  control.  It  is  quite  often  that  the  yield  stress  even 
decreaau  with  the  shear  nte,  which  further  enhaneu  the  need  for  a  low  shear  rate.  Second,  the  fluid 
Inertia  and  the  unsteady  wall  shear  stress  (which  is  higher  than  the  steady  wall  shear  stress) 
significantly  Increase  the  zero-field  damping  force.  A  flow-mode  damper  genetmly  hu  higher  shear 
nte  and  icoelendon,  eipedilly  with  a  large  A^Abr,  than  a  shear-mode  damper. 

In  the  analvsis,  the  cylinder  component  of  a  damper  Is  fixed  while  the  piston  component 
vibrates.  Both  the  components  of  many  practical  dampers  can  move.  In  most  of  these  cases,  there  is 
a  lirge  disptriw  in  the  amplitudes  of  the  motion  for  the  two  components,  and  the  analysis  in  this 
study  is  valid.  Otherwise,  more  appropriate  boundary  conditions  have  to  be  used  for  a  finer  analysis 
although  the  general  conclusion  from  this  study  may  still  apply,  at  leut  to  the  first  order  accuracy. 
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For  the  purpose  of  generality,  the  assigned  vibration  in  this  study  is  sinusoidal.  Road  surface 
profiles  for  an  automobile  damper,  for  example,  is  not  so.  Through  the  use  of  the  power  spectral 
density  function,  however,  one  can  present  the  surface  profiles  in  forms  of  sine  wave. 


CONCLUSIONS 

In  this  study,  the  dyrumics  of  the  fluid  flow  is  included  in  the  derivation  of  the  zero<field 
damping  forces  of  the  flow*,  shear-,  and  fixed-mode  ER  daix^sers.  The  main  resulu  of  the  study  are 
u  follows: 

•  The  zero-field  pressure  drop  in  a  flow-mode  damper  is  mainly  to  overcome  the  wall  shear 

stress  at  low  nequenoies.  At  high  frequencies  (with  the  dynamic  constant  h*  greater  than 
one),  the  zero-field  preuure  counters  primarily  the  fluid  menia,  and  it  increases  with  the 
frequency  or  h*^.  It  also  has  a  significant  phue  advance  relative  to  the  vibration  at  h^h 
frequencies.  The  level  of  ER  control  decreases  with  the  redo  of  the  piston  area  to  the  £R 
valve  cross-section  ApfAn. 

•  The  zero-field  preuure  drop  in  a  shear-mode  damper  also  has  to  counter  the  fluid  inerda  at  high 

frequencies.  When  the  dynamic  constant  h*  is  greater  than  0.5,  the  zero-field  pressure 
increases  with  h*  or  the  square  root  of  the  frequency. 

•  The  behavior  of  the  zero-field  preuure  drop  in  a  mixed-mode  damper  it  between  thou  of  a 

shear-  and  a  flow-mode  damper,  and  it  is  a  funcdon  of  the  area  ndo  A^Aga.  The  mixed¬ 
mode  damper  becomes  a  shesMsode  damper  when  the  redo  Ap/An  approaches  zero. 

•  To  achieve  the  same  petfcnnanoe,  a  flow-miode  ER  damper  it  not  u  eonqpaot  and  efficient  as  a 

ahear^mode  ER  danger. 

•  VHith  the  same  ER  dancing  force,  a  mixed-mode  damper  ia  more  compact  than  a  shear-mode 

damper.  However,  the  BMod-aaode  dan^  does  not  have  at  a  low  zero-field  damping  force 
at  the  ahea^mode  damper. 
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APPENDIX  B: 

Drawings  of  the  ER  Damper  for  the  HMMWV 
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APPENDIX  C: 

The  Quarter-HMMWV  Simulation  Program 


The  Simulation  Program 


C  TACSIM.FOR 
C 

C  To  simulate  the  quarter  car  model  for  TACOM  project. 

C 

C  •••••  Parameters  and  Variables 
C 

C  AS:  sprung  mass  acceleration 

C  ASA:  sprung  mass  acceleration  amplitude 

C  ASOAG:  ASA/AGA 

C  ASOAG:  ASA/AGA  array 

C  AU:  unsprung  mass  acceleration 

C  AUA;  unsprung  mass  acceleration  amplitude 

C  AUOAG:  AUA/AGA 

C  AUOAGA:  AUA/AGA  array 

C  AG:  ground  acceleration 

C  AGA:  ground  acceleration  amplitude 

C  ANGLEAi  Screw-nut  lead  angle  (e.g.  45  deg.) 

C  ANGTHR:  Thread  angle  for  Acme  threads  CW.5  deg) 

C  DEG:  phase  angle  in  degree 

C  DLS:  (dynamic)  spring  compression 

C  DLSS:  static  spring  compression 

C  DLT:  (dynamic)  tire  spring  compression 

C  DLTMAXi  maximum  tire  spring  compression 

C  DLTS:  static  tire  spring  compression 

C  DPSCRE:  Screw-nut  pitch  diameter 

C  DT:  time  step  size 

C  DTHETA:  phase  angle  step  size 

C  ERH:  ER  valve  gap  hight 

C  ERRHO:  ER  fluid  density 

C  ERRI;  ER  plate  inner  radius 

C  ERRO:  ER  plate  outer  radius 

C  ETA:  plastic  viscosity 

C  F,  FF:  Road  profile  frequency;  F  single,  FF  array 
C  FDN:  Newtonian  part  of  damping  force 

C  FDER;  ER  part  of  damping  force 

C  FDERC:  Control  or  ideal  FDERC 

C  FETA:  A  factor  in  FDN  calculation,  FDN  ■  FETA*ETA*DV 
C  FIT:  Iteration  factor 

C  FS:  spring  force 

C  FSCREVH:  Friction  coefficient 

C  FST:  total  shock  force 

C  FT:  F.t,  tire  (spring)  force 

C  FTY;  A  factor  in  FDER  calculation,  FDER  -  FTY*TY 

C  G:  Gravitational  acceleration 

C  13:  Output  interval 

C  ICONTR:  Control  strategy: 

C  0  --  No  ER  control. 


C-1 


(input) 

(in-data) 


(input) 

(input) 


(input) 

(input) 

(input) 

(input) 

(input) 


(input) 


(input) 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


1- -  Skyhook 

2- -  Conventional 
ICYCLE;  cycle  counter 
ICYCLM:  maximum  for  ICYCLE 

IMAX:  CIMAX-1)  no.  of  steps  in  each  cycle 

ITAU  -  1  if  DT/TAUER  >1.0 
-  0  otherwise 
NPLT:  No.  of  plates 

OMEGA:  angular  frequency 

RFT:  Ratio  of  the  translational  damping  force  to  torque 

RK:  k  in  equation  CW3) 

RKS:  spring  constant. 

RKSC:  k_c,  compression  constant  in  eq  w4 

RKST:  k_t,  tension  constant  in  eq  w4 

RKT:  tire  spring  constant 

RLS:  (dynamic)  spring  length 

RLSS:  static  spring  length 

RMS:  sprung  mass 

RMU:  unsprung  mass 

RTR:  Ratio  of  the  translational  velocity  to  rotational 

or  angular  velocity 
T:  time 

TAUER:  ER  fluid  time  constant 
THETA:  phase  angle 

TY:  The  Bingham  stress  in  TY  -  TYC*CVC)atyN 

TYCi  A  constant  in  TY  -  TYC*CVC)atyn 
TYLIM:  TY  at  VCLIM 

TYN:  The  exponent  in  TY  -  TYC*CVC)atyn 

VC:  control  voltage 

VCLIM  the  limit  of  VC 

VG:  gournd  velocity 

VGA:  ground  velocity  amplitude 

VS:  sprung  velocity 

VU:  unsprung  velocity 

VV  -VSCI)*VSU 

XG:  ground  displacement 

XGA;  ground  displacement  amplitude 

XS:  sprung  mass  displacement 

XU:  unsprung  mass  displacement 

IMPLICIT  REAL*8CA-H,0-Z) 

DIMENSION  XSC5000) ,XUC5000) ,VS(5000),VUC5000) 

DIMENSION  ASC5000),AUC5000> 

DIMENSION  FSC5000) , FDERC5000) , FDNC50005 , DEGC5000) 
DIMENSION  FSTC5000),FTC5000) 

DIMENSION  FF(50)  ,ASOAGAC.50)  ,AUOAGAC50) 

OPEN  CS.FILE-'tacsim.S') 

OPEN  C6,FILE-'tacsim.6') 

READC5,*)  ANGLEA 
READC5,*)  DLTMAX 
READC5,*)  DPSCRE 
READC5,*)  ERH 
READC5,*)  ERRHO 
READC5,*)  ERRI 
READCS,*)  ERRO 


(input) 

(input) 


(input) 


(input) 

(input) 

(input) 

(input) 

(input) 


(input) 

(input) 


(input) 

(input) 

(input) 

(input) 


(input) 


C-2 


READ(5,*)  ETA 
READC5,*)  FIT 
READ(5,*)  FSCREW 
READ(5,*)  13 
READC5,*)  ICONTR 
READC5,*)  ICYCLM 
READ(5,*)  IMAX 
READ(5,*)  NPLT 
READ(5,*)  RK 
READ(5,*)  RKS 
READ(5,*)  RKSC 
READC5,*)  RKST 
READ(5/)  RKT 
READ(5,*)  RMS  - 
READ(5,*)  RMU 
READC5,*)  TAUER 
READC5,*)  TYC 
READ(5,*)  TYN 
READ(5,*)  VCLIM 
READC5,*)  XGA 
READ(5,*)  NF 
DO  5  IFF=1,NF 
READC5,*)  FF(IFF) 

5  CONTINUE 

‘  WRITEC6,*)  'ANGLEA=' .ANGLEA 
WRITE(6 , *)  ' DLTMAX= ' , DLTMAX 
WRITEC6 , *)  ' DPSCRE= ’ , DPSCRE 
WRITE(6,*)  'ERH  =’,ERH 
WRITEC6,*)  'ERRHO  =',ERRHO 
WRITECe,*)  'ERRI  =',ERRI 
WRITEC6,*)  'ERRO  =',ERRO 
WRITE(6,*)  'ETA  =',ETA 
WRITE(6,*)  'FIT  =',FIT 
WRITEC6 , ' FSCREW= ' , FSCREW 
WRITE(6,*)  '13  =',I3 

WRITECe , *)  ' ICONTR= ' , ICONTR 
WRITEC6 , *)  ' ICYCLM= ' , ICYCLM 
WRITE(6,*)  'IMAX  =' ,IMAX 
WRITEC6,*)  'NPLT  =' ,NPLT 
WRITE(6,*)  'RK  =',RK 
WRITE(6,*)  'RKS  =',RKS 
WRITE(6,*)  'RKSC  =' ,RKSC 
WRITEC6,*)  'RKST  =' ,RKST 
WRITE(6,*)  'RKT  =' ,RKT 
WRITE(6,*)  'RMS  =' ,RMS 
WRITEC6,*)  'RMU  =' ,RMU 
WRITEC6,*)  'TAUER  =', TAUER 
WRITE(6,*)  'TYC  =' ,TYC 
WRITEC6/)  'TYN  =' ,TYN 
WRITE(6,*)  'VCLIM  =' , VCLIM 
WRITE(6,*)  'XGA  =' ,XGA 
WRITECe,*)  'NF  =',NF 
DO  7  IFF=1,NF 

7  WRITE(6,*)  'F'.IFF,'  =  ',FF(IFF) 


Q  Initialization 


WRITE(6,*)  '  ' 

WRITEC6,*)  'Some  output  ' 

G-9. 89665 
PI-3. 14159 

C  thread  angle  14.5  deg. 

ANGTHR-14.5*PI/180.0 
C  lead  angle  unit  conversion 
ANGLEA-ANGLEA*PI/189.0 
C  Screw-nut  conversion  ratios 

RFT-CDPSCRE/2.0)*CDrOS(ANGTHR)*DTANCANGLEA)- 
+  FSCREW)/CFSCREW*DTANCANGLEA)+DCOSCANGTHR)) 

RTR-1 .0/CCDPSCRE/2 . 0)'*DTANCANGLEA)) 

WRITEC6,*)  'RFT  «',RFT 
WRITEC6,*)  'RTR  -',RTR 

C  FTY  defined  in  FDER  «  FTY*TY 

FTY-C4.*PI*NPLT/3.)*CERR0**3-ERRI**3)/RFT 
writeCG/)  'FTY  -  FTY 

C  FETA  defined  in  FDN  -  FETA*CDV*ETA) 

FETA»CPI*NPLT*RTR/ERH)*CERR0**4-ERRI**4)/RFT 
writeCe,*)  'FETA  -',FETA 

TYLIM-TYC'»VCLIM*'*TYN 
WRITECG,-*)  'TYlim  -'.TYLIM 

C  Spring  stuff 

DLSS-RMS»G/RKS 

DLTS-CRMS+RMU)*G/RKT 

WRITEC6,'*)  'static  spring  deflection  -'.DLSS 
WRITE(6,*)  'static  tire  deflection  -',DLTS 

^  pf^EQUENCY  ITERATION 

DO  800  IFF-1,NF 

F-FFCIFF) 

WRITEC6,*)  ■  ' 

WRITECe,*)  '  ****••••*•***•*  F-  '.F 

C  Angular  frequency 
OMEGA-2. 0'*PI*F 

C  Ground  acceleration  amplitude 
AGA-0MEGA**2*XGA 

C  Ground  velocity  amplitude 

VGA-OMEGA*XGA 
WRITECG,-*)  'AGA  ,AGA 

WRITE(6,*)  'VGA  ,VGA 

C  phase  angle  step  size  (in  Radium) 

DTHETA-2.0*PI/CIMAX-1.0) 

C  time  step  size 

DT-DTHETA/OMEGA 
WRITEC6,*)  'DT  DT 

CONST-DT/TAUER 

WRITEC6,*)  'DT/TAUER  -' .CONST 
ir(CONST.GT.1.0)  THEN 
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ITAU^l 
ELSE 
ITAl'-e 
END  IF 

(;«**««**«*•*•**«*«*  Marching  in  Time  Domain 
I-l 

DEGCl)-0.e 

XUCI)-t).0 

XSCI)-0.0 

VSCI)-0.0 

VUCI)-0.0 

AS(I)-0.0 

AU(I>0.0 

FSCI)-0.0 

FT(I)-0.0 

FDNCI)-0.0 

FDERCI>0.0 

FSTCI)-0.0 

TYOI.,D-0.0 

ICYCLE-0 
ASA-9999.999 
AUA-999g.999 
40  ICYCLE-ICYCLE+1 


IFCICYCLE.GT.l)  THEN 

DASA-DABSC(ASA-ASAOLD)/ASACLD) 

I-l 

XUCI>XUCIMAX) 

XSCI>XS(IMAX) 

VSCI)-VSCIMAX) 

VUCI)-VU(IMAX) 

ASCI)-AS(IMAX) 

AUCD-AUCIMAX) 

FSCD-FSCIMAX) 

FTCI)-FT(IMAX) 

FDNCD-FDNCIMAX) 

FDER(I)-FDERCIMAX) 

FSTCD-FSTCIMAX) 

I FCDA5A.lt. 0.0001)  GO  TO  100 
END  IF 

ir(ICYCLE.GT.ICYCLM)  GO  TO  100 

ASAOLD-ASA 

AUAOLD-AUA 

ASA-i0.0 

AUA-0.0 

VCMAX-0.0 

TYMAX-0.0 

IMAXMl-IMAX-1 

C  NUMBER  OF  ITERATIONS 
NIT-0 

C  AVERAGE  CONVERGENCE 
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ACIT-0.0 
DO  90  I-2,IMAX 

theta-ci-i.0)*dtheta 

DEGCD-THETA*  180. 0/PI 


C  ground  Motion 

XG-XGA*DSINCTHETAD 
C  ground  velocity 

VG-CMEGA*XGA*DCOSCTHETA) 

xsci>xsci-i) 

XUCI)-XUCI-1)  ■ 

VSCI>VSCI-1) 

VUCD-VUCI-I) 

ASCD-ASCI-l) 

AUCI>AUCI-1) 


Q  *«*•«•**««««*****•«*  iteration 

DO  70  IT-1,20 

C  spring  force 

DLS-DLSS+XU(I)-XSCI) 

FSCI>RKS*DLS 
C  relative  velocity 
OV-VUCD-VSCD 

C  •••**  Newtonian  damping  force  CFDN) 

C  writeCe,*)  'stop  1' 

FDNCI)-FETA*ETA*DV 

C  '*''"**  ER  damping  force  (FDER) 

C  --  No  ER  control 

IFCICONTR.EQ.0)  THEN 
FDERCI>0.0 
END  IF 

C  --  Control  Strategy  1 
C  Skyhook 

IFCICONTR.EQ.l)  THEN 
VSU-VSCI)-VUCI) 

VV-VSCI)*VSU 
IF(VV.GT.0.0)  THEN 

C  FDER  and  VS  have  the  opposite  sign. 

C  Control  voltage 

VC-CDABSCRK*VSCI)))**(1 . /2 . 5) 

IFCVC.GT.VCLIM)  VC-VCLIM 
IFCITAU.EQ.0)  THEN 

TY-TYOLD+(DT/TAUER)*CTYC'*VC'**TYN-TYOLD) 

ELSE 

TY-TYC*VC**TYN 
END  IF 

fderci)-cvsci)/cdabscvsci))+0.0000001))*fty*ty 
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ELSE 

FDERCI)-8.0 
END  IF 
END  IF 

C  “  Control  Strategy  2 

C  Conventional  Damper 

IFCICONTR.EQ.2)  THEN 
VSU-VSCI)-VUCI) 

I FCVSU.lt. 0.0)  THEN 

C  COMPRESSION  (damping  force  on  the  sp^'ung  mass  pointing  up,  FDER  >  0) 
VC-CDABS(RKSC*VSU))* •(! . /2 . 5) 

IF(VC.GT.VCLIM)  VC-VCLIM 
IF(ITAU.EQ.0)  THEN 

TY«TYOLD+CDT/TAUER)*CTYC*VC**TYN-TYOLD) 

ELSE 

TY-TYC*VC**TYM 
END  IF 

FDERCI)— CVSU/(DABSCVSU)+0.0000001))*FTY*TY 
ELSE 

C  TENSION  (damping  force  on  the  sprung  mass  pointing  down,  FDER  <  0) 
VC-(DABS(RKST*VSU))**(1 . /2 . 5) 

IF(VC.GT.VCLIM)  VC-VCLIM 
IF(ITAU.EQ.0)  THEN 

TY-TYOLD+CDT/TAUER)*(TYC*VC**TYN-TYOLD) 

ELSE 

TY-TYC*VC**TYN 
END  IF 

FDER(I)— (VSU/(DABS(VSU)+0 . 0000001))*  FTY*TY 
END  IF 
END  IF 

C  Control  electric  field  strength 
C  Total  strut  force 

FSTCI)-FDER(I)+FDN(I)+FS(I) 


C  ••••*•**••  Integration 

C  Unsprung  mass  dynamic.s 

C  tire  compression  (positive  when  compressed) 

DLT-DLTS-XU(I)+XG 
IF(DLT.GT.DLTMAX)  THEN 

WRITE(6,*)  'STOP,  TIRE  BOHOMED  OUT.  ICYCLE,  I-' ,ICYCLE,  ’  ,  ’  ,I 
CLOSECS) 

CLOSEC6) 

STOP 
END  IF 

IF(DLT,GT.0.0)  THEN 
FT(I)-RKT*DLT 
ELSE 

FT(I)-0.0 
END  IF 

C  unsprung  mass  acceleration 
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CONST-C-FSTCI)+FTCI))/RMU-G 

auci>ci.0-fit)»au(i)+fit*const 

C  velocity  and  displacement 

VUCI>VUCI-1)+0.5*CAUCI)+AUCI-1))*DT 
XUCI>XUCI-l)+0 . 5*CVUCI)+VU(I-1))*DT 

C  ••••  Sprung  mass  dynamics 

C  sprung  mass  acceleration 
CONST-FSTCI)/RMS-G 

CIT-DABSCCONST-ASCI))/(DABSCASCI))+0. 0000000001) 
ASCD-Cl .  0-  FIT)*ASCI)+FIT*C0NST 

C  velocity  and  displacement 

VSCI)-VSCI-1)+0.5*CAS(I)+ASCI-1))*DT 

XSCI)-XSa-l)+0.5*CVSCI)+VSCI-l))*DT 

IFCn.GT.2)  THEN 

I FCCIT.lt. 0.001)  GO  TO  71 
END  IF 

70  CONTINUE 

71  CONTINUE 

NIT-NIT+IT-1 

ACIT-ACIT+CIT 

C  absolute  value 

AAS-DABS(ASCI)) 

IFCAAS.GT.ASA)  THEN 
ASA-AAS 
END  IF 

AAU-DABSCAUCD) 

IFCAAU.GT.AUA)  THEN 
AUA-AAU 
END  IF 


C  Update 


IFCVC.GT.VCMAX)  VCMAX-VC 
IFCTY.GT.TYMAX)  TYMAX-TY 

TYOLD-TY 

90  CONTINUE 

GO  TO  40 

100  CONTINUE 

ASOAG-ASA/AGA 
AUOAG-AUA/AGA 
WRITEC6,*)  •' 

WRITECe,*)  'ASA/AGA  -',ASOAG 
WRITE(6,*)  'AUA/AGA  ,AU0AG 
ICYCLE~ICYCLE-1 
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WRITEC6,*)  'ICYCLE  -'.ICYCLE 
WRITEC6,*)  'VCMAX  ,VCMAX 

WRITEC6,*)  'TYMAX  JYMAX 

ASOAGACIFF)-ASOAG 

AU0AGAC1FF)-AU0AG 

NIT-NIT/IMAXMl 

WR1TEC6,*)  'AVERAGE  NUMBER  OF  ITERATION  -  ',NIT 
ACIT-ACIT/IMAXMl 

WRITEC6,*)  'AVERAGE  CONVERGNECE  (%)  -  ' ,ACIT 

600  FORMATCE12.5) 

IFCNF.EQ.l)  THEN 

WRITEC6,'»)  •  ' 

WRITEC6,'»)  'DEG  -' 

DO  120  I-1,IMAX,I3 
120  WRITEC6,600)  DEGCD 

WRITEC6,'*)  ’  ' 

WRITE(6,'»)  'XS 
DO  122  I«1,IMAX,I3 

122  WRITEC6,600)  XSCD 

WRITEC6,*)  '  ' 

WRITEC6,*)  'XU 
DO  123  I-1,IMAX,I3 

123  WRITEC6,600)  XUCI) 

WRITEC6,'*)  '  ' 

WRITEC6,*)  'VS  -• 

DO  124  I-1,IMAX,I3 

124  WRITEC6,600)  VSCI) 

WRITECe/)  '  ' 

WRITECe,'*)  'VU 
DO  125  I-1,IMAX,I3 

125  WRITEC6,600)  VUCD 

WRITEC6,'*)  '  ' 

WRITEC6,*)  'AS 
DO  322  I-1,IMAX,I3 
322  WRITEC6,600)  AS(I) 

WRITEC6,*)  '  ' 

WRITEC6,*)  'AU  -' 

DO  324  I-1,IMAX,I3 
324  WRITEC6,600)  AUCI) 

WRITEC6,*)  '  ' 

WRIT£(6,*)  'FS 
DO  126  I-1,IMAX,I3 

126  WR1TEC6,600)  FS(I) 

WRITECe,')  '  ' 

WRITE(6,*)  'FT  -' 
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DO  127  I-1,IMAX,I3 

127  WRITE(6,600)  FTCD 

WRITE(6,*)  '  ' 
WRITEC6,*)  'FDN 
DO  128  I-1,IMAX,I3 

128  WRITEC6,600)  FDNCI) 

WRITEC6,*)  '  ' 
WRITEC6,*)  'FDER 
DO  130  I-1,IMAX,13 
130  WRITEC6,600)  FDERCI) 

WRITEC6,*)  '  '  ■ 
WRITEC6,*)  'FST 
DO  136  I»1,IMAX,I3 
136  WRnEC6,60C)  FSTCI) 

END  IP 

800  CONTINUE 

WRITEC6,*)  '  ' 
INRITEC6,*)  'ASA/AGA' 

DO  810  IFF-1,NF 
810  WRITECe,*)  ASOAGACIFF) 

WRITEC6,*)  '  ' 
WRITEC6,*)  'AUA/AGA' 

DO  820  IFF-l.NF 
820  WRITECe,*)  AUOAGA(IFF) 

CL0SEC5) 

CL0SEC6) 

STOP 

END 


An  Input  Data  File 

17.6568, 

ANGLEA: 

Screw-nut  lead  angle  (e.g.  45  deg.) 

0.25, 

DLTMAX: 

maximum  tire  spring  compression 

(input) 

0.0254, 

DPSCRE: 

Screw-nut  pitch  diameter,  1"  pitch,  1" 
lead 

(input) 

0.001, 

ERH: 

ER  valve  gap  hight 

(input) 

1000.0, 

ERRHO: 

ER  fluid  density 

(input) 

0.0285, 

ERRI: 

ER  plate  inner  radius,  Di  -  2.244",  Ri  - 
.0285  m 

(input) 

0.0483489, 

ERRO: 

ER  plate  outer  radius,  Do-3.807", 

Ro-0. 04835  m 

(input) 

0.222, 

ETA: 

plastic  viscosity,  100cP  -  0.1  N.s/m^2 
(input) 

0.2, 

FIT 

iteration  factor,  - 

0.033, 

FSCREW: 

Friction  coef  -  0.033  For  a  efficiency  of 

89*  with  1"X1" 

5, 

13. 

Output  interval 

(input) 

1, 

ICONTR: 

Control  strategy  - 

(input) 

C-IO 


10, 

ICYCLEM: 

Maximum  cycle  number 

361, 

IMAX: 

(IMAX-1)  no.  of  steps  in  each  cycle 

40, 

NPLT: 

No.  Of  plates 

590550000.0, 

RK: 

k  in  equation  w3 

163300.0, 

RKS: 

spring  constant 

315000000.0, 

RKSC! 

k_c,  compression  constant  in  eq  w4 

630000000.0, 

RKST: 

k_t,  tension  constant  in  eq  vi/4 

463800.0, 

RKT: 

tire  spring  constant 

998.0, 

RMS: 

sprung  mass  2200  lb  -  998  kg 

181.0, 

RM 

U:  unsprung  mass  400  Ibm  >•  181  kg 

0.0002, 

TAUER 

0.000095, 

TYC 

2.0, 

TYN 

3000.0 

VCLIM  - 

LIMIT  FOR  VC  -  3000  V. 

0.01, 

XGA: 

ground  displacement  amplitude. 

25, 

NF 

le-1. 

FI 

2e-l, 

F2 

Se-l, 

8e-l» 

le-f0) 

1.5e-»-0t 

1.86+0, 

1.9e+0, 

2e+0, 

2.5e+0, 

3e+0, 

4e+0, 

6e+0, 

7«+0, 

8e+0, 

8 ■ 56+0 , 

96+0, 

9.56+0, 

l6+l, 

1.2e+l, 

1.4e+l, 

1.6e+l, 

1.8e+i, 

2e+l, 

2.5e+l, 


(input) 

(input) 

(input) 

(input) 

(input) 

(input) 

(input) 

(input) 

(input) 


(input) 
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APPENDIX  D: 

Data  Plots  from  Tests  of  the  ER  Damper 
on  the  Quarter-HMMWV  Ride  Simulator 


This  appendix  presents  a  complete  set  of  data  plots  from  the  testing  of  the  ER  damper 
on  the  quarter-HMMWV  ride  simulator. 

Figure  D-1  shows  the  modified  road  profiles  used  to  drive  the  simulator.  This  profiles 
were  derived  from  the  proving  grounds  road  profiles  which  were  provide  to  UMTRI  by 
TACOM.  The  modifications  were  essentially  the  necessary  high-pass  filtering  required  to 
reduce  the  vertical  elevation  of  the  road  profiles  to  within  the  capability  of  the  simulator. 

Figure  D-2  through  D-2S  present  results  of  testing  displayed  as  transfer  function  gains 
in  the  frequency  domain.  Figures  D-2  through  D-9  show  the  gains  of  sprung  mass 
acceleration,  relative  to  road  surface  acceleration,  as  observed  on  the  eight  road  surfaces, 
respectively.  Each  figure  contains  three  plots,  one  each  for  the  gains  under  three  operating 
conditions.  These  were  (1)  HR  damper  installed  with  control  voltage  off,  (2)  ER  damper 
installed  and  the  skyhook  controller  in  use,  and  (3)  conventional  HMMWV  shock  absorber 
installed. 

Figures  D-10  through  D-1 7  is  display  comparing  the  sprung  mass  acceleration  gains 
observed  using  the  conventional  shock  absorber  with  the  same  gain  observed  when  the 
ER  damper  was  controlled  to  simulate  a  conventional  damper.  Again,  there  are  eight  figures 
associated  with  the  eight  roads,  respectively. 

Finally,  figures  D-18  through  D-25  show  the  tire  forces  gains,  that  is:  (tire  force)/(ruad 
surface  acceleration).  The  eight  figures  are  for  the  eight  roads,  respectively.  The  three  plots 
are  for  the  same  three  operating  conditions  described  for  figures  D-2  through  D-17. 
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Figure  D.l.  Road  profiles  used  In  quarter-HMMWV  testing 
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Figure  D.l  (continued).  Roud  proflies  used  in  quarter-HMMWV  testing 
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Gain:  AzSPRUNQ  /  AzROAD  K  -  0  (off) 
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Figure  D.2.  Sprung  mass  acceleration  transfer  function  gains  for  APG  09 
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Figure  D.3.  Sprung  mass  (acceleration  transfer  function  gains  for  APG  29 
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Gain:  AzSPRUNG  /  AzBOAD  K  «  0  (o«) 
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Figure  D.4.  Sprung  mass  aeceleration  transfer  function  gains  for  APG  37 
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Figure  D.5.  Sprung  mass  acceleration  transfer  function  gains  for  CHV  01 
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Gain:  AzSPRUNG  /  AzROAD  Conventional  Army  Shock 
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Figure  D.6.  Sprung  mas.s  acceleration  transfer  function  gains  for  CHV  06 
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Figure  D.7.  Sprung  mass  acceleration  transfer  function  gains  for  FTKN  31 
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Figure  D.8.  Sprung  mass  acceleration  transfer  function  gains  for  FTKN  34 
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Figure  D.9.  Sprung  muss  accciet  jUon  transfer  ftinction  gains  for  LET  5 


Qain;  AzSPRUNG  /  AzROAD  Simulatad  Convantlonal  (Kb8«6/16«6) 


Qaln;  AzSPRUNG  /  AzROAD  ConvtntiQnal  Army  Shook 


Jl 

u 

m 

_ 

■ 

m 

■HS 

iSfSBSRSKII 

0  5  10  IB 

Frequency  •  Hz 


Figure  D.IO.  Sprung  mass  acceleration  transfer  function  gains  measured 
with  simulated  and  actual  conventional  shock  absorber  for  APG  09 
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Figure  D.ll.  Sprung  mass  acceleration  transfer  function  gains  measured 
with  simulated  and  actual  conventional  shock  absorber  for  APG  29 
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Figure  D.12.  Sprung  mass  acceleration  transfer  function  gains  measured 
ivith  simulated  and  actual  conventional  shock  absorber  for  APG  37 
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Figure  D.13.  Sprung  mass  acceleration  transfer  function  gains  measured 
with  simulated  and  actual  conventional  shocl<  absorber  for  CHV  01 
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Figure  D.14.  Sprung  mass  acceleration  transfer  ftinctiun  gains  measured 
with  simulated  and  actual  conventional  shock  absorber  for  CHV  06 
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Figure  D.i6.  Sprung  mass  acceleration  transfer  ftinctton  gains  measured 
with  simulated  and  actuai  conventional  shuck  absorber  for  FTKN  34 
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Figure  D.17.  Sprung  mass  acceleration  transfer  function  gains  measured 
with  simulated  and  actual  conventional  shock  absorber  for  LET  5 
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Fisuru  D.18.  Tire  force  transfer  function  gains  for  APG 
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Figure  D.20.  Tire  force  transfer  function  gains  for  APG  37 
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Figure  D.21.  Tire  force  transfer  function  gains  for  CHV  01 
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Figure  D.22.  Tire  force  transfer  function  gains  for  CHV  06 
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Figure  D.23.  Tire  force  transfer  function  gains  for  FTKN  31 
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Figure  D.24.  Tire  force  transfer  function  gains  for  FTKN  34 
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Figure  D.25.  Tire  force  transfer  function  gains  for  LET  5 


